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1  INTRODUCTION 


JASON  has  been  tasked  by  the  Army  to  review  recent  progress  in  in¬ 
frared  detector  technology  and  to  update  the  earlier  JASON  report  “Imaging 
Infrared  Detectors”  (JSR-91-600).  To  this  end,  a  workshop  was  conducted 
during  the  1997  Summer  Study  with  distinguished  participants  drawn  from 
industry,  academia,  and  government,  listed  in  Appendix  D.  Through  dis¬ 
cussions  with  workshop  participants  and  the  JASONs  an  updated  overview 
of  the  state  of  imaging  infrared  detector  technology  was  developed,  and  is 
presented  in  the  body  of  this  report  as  follows:  Section  2  (Image  Character¬ 
ization);  Section  3  (Uncooled  Microbolometer  Arrays);  Section  4  (Quantum 
Well  Infrared  Photodetectors);  Section  5  (HgCdTe  Photodetector  Arrays); 
and  Section  6  (New  Developments).  Appendices  to  the  original  report  con¬ 
tained  primers  on  infrared  basics,  which  have  been  proven  to  be  useful.  The 
current  report  contains  updated  versions:  Appendix  A  (Fundamentals  of  IR 
Detector  Performance);  Appendix  B  (Imaging  Bolometer  Arrays);  and  a  new 
addition,  Appendix  C  (Thermoelectric  Materials). 

Imaging  infrared  detector  technology  has  advanced  rapidly.  Important 
developments: 

•  Driven  by  the  need  for  image  recognition,  large  staring  long  wave  in¬ 
frared  (LWIR)  arrays  have  been  developed:  320x240  pixel  uncooled 
bolometer,  and  cooled  quantum  well  infrared  photodetector  (QWIP), 
and  cooled  HgCdTe  detector  arrays  are  available  as  deliverable  prod¬ 
ucts;  640x480  pixel  arrays  have  been  demonstrated  for  the  QWIP 
and  HgCdTe  technologies.  Larger  arrays  up  to  high  density  television 
(HDTV)  are  planned  for  all  three  approaches. 
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•  Uncooled  bolometer  and  cooled  QWIP  cameras  with  excellent  perfor¬ 
mance  are  commercially  available,  both  inside  the  U.S.  and  abroad. 

•  Two  “color”  (e.g.  LWIR  and  MWIR)  detector  arrays  are  under  develop¬ 
ment,  with  two  color  pixel  demonstrations  for  the  QWIP  and  HgCdTe 
technologies. 

•  Continued  improvements  in  Si  readout  integrated  circuits  (ROIC)  have 
led  to  better  uniformity  and  sensitivity. 

•  Methods  for  mechanical  strain  relief  between  the  detector  array  and  the 
Si  ROIC  have  been  developed  -  e.g.  using  thinned  detector  materials 
and  ROIC  strain  matching  -  for  improved  reliability. 

Recent  advances  in  infrared  camera  technology  will  have  substantial 
impact  on  the  Army: 

•  Uncooled  IR  cameras  now  offer  performance  traditionally  associated 
with  cooled  HgCdTe  units,  but  at  relatively  low  prices.  They  are  light 
weight,  have  low  power  consumption,  and  produce  excellent  images 
with  good  sensitivity  (noise  effective  temperature  difference  NETD  ~ 
70  mK  with  f/ 1  optics) .  In  the  military  uncooled  IR  cameras  will  be  the 
natural  choice  for  small  unit  operations,  surveillance,  and  even  some 
targeting  applications.  These  cameras  have  many  possible  commercial 
applications  -  such  as  security,  fog  penetration,  and  night  driving  - 
which  may  drive  future  development  work.  The  international  avail¬ 
ability  of  uncooled  technology  will  soon  change  “owning  the  night”  to 
“sharing  the  night.” 

•  QWIP  technology  has  been  used  to  make  small,  portable  LWIR  cam¬ 
eras  using  miniaturized  coolers.  Advantages  of  QWIP  arrays  are  ex¬ 
cellent  image  quality  due  to  high  uniformity,  and  a  mature  material 
technology  permitting  the  production  of  large  arrays  (640x480  pixels) 
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and  complex  pixel  structures  for  multi-color  and  polarization  sensitive 
detectors.  However,  their  quantum  efficiency  is  limited  to  relatively 
small  values  by  fundamental  considerations,  requiring  cooling  to  70  K 
and  below  for  best  sensitivity. 

•  Large  staring  HgCdTe  LWIR  cameras  with  640x480  pixels  have  been 
demonstrated.  These  instruments  have  excellent  image  quality  and 
high  sensitivity,  and  would  be  the  natural  choice  for  sensitive  IR  tele¬ 
scopes  to  be  used  for  surveillance  and  targeting  at  long  ranges.  These 
arrays  are  impressive  demonstrations  of  progress  in  HgCdTe  technol¬ 
ogy.  But  they  are  expensive,  and  their  use  will  be  limited  to  the  most 
demanding  applications  while  relatively  inexpensive  uncooled  cameras 
will  suffice  for  most  applications.  The  fabrication  technology  for  LWIR 
HgCdTe  detectors  is  still  not  fully  mature,  and  difficulties  remain. 
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2  IMAGE  CHARACTERIZATION 


The  requirements  placed  on  infrared  cameras  by  Army  applications  can 
be  quite  different  than  airborne  or  strategic  applications.  The  central  func¬ 
tion  of  infrared  cameras  for  Army  use  is  the  detection  and  recognition  of 
desired  images  -  soldiers,  fixed  and  mobile  military  equipment  -  in  the  pres¬ 
ence  of  spatial  clutter  -  grass,  trees,  rocks,  civilian  cars  and  trucks.  Camera 
noise  and  nonuniformity  also  produce  a  type  of  clutter.  The  detection  of 
small,  featureless  targets,  such  as  missile  exhausts,  is  also  important,  but 
does  not  play  the  dominant  role  it  does  in  space.  For  the  separation  of  im¬ 
ages  from  clutter,  the  spatial  characteristics  of  the  image  and  the  detector 
array  are  as  important  as  the  temporal  characteristics. 

The  push  toward  larger  detector  arrays  with  more  pixels  in  recent  years 
has  been  driven  by  the  realization  that  many  pixels  on  target  are  needed  to 
separate  important  images  from  clutter.  A  single  bright  pixel  is  probably 
not  meaningful  in  the  presence  of  clutter,  unless  it  has  additional  spatial 
characteristics  (e.g.  rapid  motion)  to  identify  it.  An  image  of  a  military 
target,  such  as  a  tank,  with  enough  pixels  for  recognition,  can  be  separated 
from  the  background  without  additional  information.  Even  more  pixels  are 
desirable  to  identify  the  target,  distinguish  friend  from  foe. 

For  terrestrial  use,  balanced  figures  of  merit  which  give  equal  weight  to 
spatial  and  temporal  “noise,”  are  needed  to  evaluate  infrared  imager  per¬ 
formance  .  Traditionally,  infrared  detectors  have  been  evaluated  primar¬ 
ily  by  their  temporal  noise  through  figures  of  merit  like  the  noise  effective 
power  (NEP)  and  the  NETD  of  a  single  pixel.  Spatial  “noise”  such  as  non¬ 
uniformity  and  fixed  pattern  noise  has  not  been  given  equal  emphasis,  even 
though  spatial  “noise”  can  dominate  the  subjective  image  quality.  Better, 
balanced  figures  of  merit  are  needed  to  evaluate  and  compare  the  perfor¬ 
mance  of  infrared  imagers.  A  possible  example  is  the  spatio-temporal  power 
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spectrum  S(k,u )  where  k  is  the  two-dimensional  spatial  wavevector  across 
the  image  plane,  and  u  is  the  temporal  frequency.  The  dependence  of  the 
power  spectrum  S(k,  u)  on  the  magnitude  of  k  could  distinguish  white,  power 
law  (l//a),  and  periodic  spatial  noise,  and  plots  of  constant  S(k,  oj)  surfaces 
in  the  three  dimensional  k,u  space  would  uncover  and  characterize  spatial 
anisotropy,  and  spatio-temporal  characteristics  like  diffusion  and  cross-talk. 
The  use  of  spatio-temporal  power  spectra  for  characterization  is  standard  in 
physics  and  materials  science  for  systems  with  linear  responses.  For  detec¬ 
tors  (and  for  the  field  of  computer  vision)  in  which  non-linear  effects  can  be 
important,  other  approaches  may  be  better,  and  should  be  explored.  But 
over-simple  characterizations  of  the  spatial  characteristics  of  infrared  detec¬ 
tor  arrays  (e.g.  the  “nonuniformity”  is  X%)  are  unlikely  to  be  adequate. 

Recognition  of  targets  and  threats  is  the  most  important  function  of  in¬ 
frared  cameras  for  military  users.  The  recognition  range  is  related  to  detector 
sensitivity  (NEP)  in  a  simple  manner  explained  below,  but  this  connection 
is  not  easily  summarized  in  a  single  graph.  To  understand  the  relationship, 
consider  the  problem  of  target  identification  and  recognition  for  a  tank  mov¬ 
ing  toward  the  user  from  far  away,  first  without  atmospheric  attenuation  due 
to  water  vapor,  fog  and  other  obscurants,  and  then  with  attenuation.  Figure 
1  illustrates  the  situation  for  parameters  appropriate  to  an  uncooled  LWIR 
camera.  At  long  ranges  R,  the  tank  fills  less  than  one  pixel  of  the  detec¬ 
tor  array,  the  temperature  sensitivity  increases  as  the  tank  moves  closer  as 
NETD  oc  1/R2,  and  the  target  recognition  range  Rr  increases  with  detector 
sensitivity  as  Rr  oc  1/ VNEP .  In  this  long  range  limit,  the  detection  range 
for  point-like  objects  increases  with  sensitivity,  according  to  the  conventional 
picture.  However,  once  the  tank  occupies  more  than  one  pixel  of  the  image, 
the  temperature  sensitivity  NETD  becomes  independent  of  range  R ,  as  indi¬ 
cated  in  Figure  1.  The  reason  is  that  the  1/R2  increase  in  photon  flux  from 
the  tank  at  short  ranges  R  is  cancelled  by  the  increase  in  number  of  pixels 
occupied  by  the  image.  If  one  can  see  the  tank  100  ft.  away  in  clear  weather, 
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Figure  1:  Illustration  of  the  importance  of  pixels  on  target  for  image  recog¬ 
nition,  and  of  the  behavior  of  noise  effective  temperature  difference  (NETD) 
with  range.  This  example  assumes  typical  uncooled  camera  specifications: 
100  mm  focal  length,  ~  10°  field  of  view,  and  50  pm  pixels,  which  give  ~  1 
m  spatial  resolution  at  2  km. 


one  can  also  see  and  identify  it  at  1  km,  2  km,  and  longer  ranges,  until  it 
disappears  into  a  point  in  one  pixel,  as  illustrated  in  Figure  1,  in  accord  with 
everyday  experience.  It  is  this  many-pixels-on-target  regime  that  is  most 
important  for  Army  applications.  In  this  regime  the  recognition  range  is  in¬ 
dependent  of  detector  sensitivity  (NEP)  for  signal  to  noise  ratios  S/N  i  1  in 
each  pixel;  if  a  target  can  be  seen  at  short  range,  it  can  be  seen  at  all  ranges 
until  it  disappears  into  a  point.  Detector  sensitivity  becomes  important  in 
setting  the  intensity  threshold  for  recognition  at  any  range.  Less  infrared 
flux  is  received  from  the  target  for  small  pixel  sizes,  such  as  the  25x25  pm2 
pixels  now  entering  production  to  increase  array  size  within  chip  area  limits, 
and  for  smaller  numerical  aperture  optics,  required  by  size  and  weight  con¬ 
siderations  for  infrared  telescopes.  Thus  the  effect  of  detector  noise  for  clear 
atmospheric  conditions  is  to  place  lower  bounds  on  pixel  size  and  numerical 
aperture,  not  to  limit  recognition  range. 
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Often  targets  are  obscured  by  absorption  and  scattering  due  to  water 
vapor,  fog,  smoke,  and  other  aerosols,  particularly  under  battle  conditions. 
Under  these  circumstances  the  transmitted  infrared  photon  flux  falls  more 
rapidly,  and  can  be  often  characterized  by  an  exponential  attenuation  length 
Ra-  The  attenuation  length  varies  greatly  with  atmospheric  conditions.  To 
illustrate  the  effects  on  recognition,  again  consider  a  tank  approaching  from 
far  away,  but  under  foggy  conditions.  Initially  its  image  size  in  the  camera 
is  less  than  one  pixel,  and  the  temperature  sensitivity  falls  off  with  range 
as  NEPD  oc  exp(— R/RA)/R2.  When  the  image  size  reaches  one  pixel  (at 
range  R\  determined  by  camera  optics),  the  tank  may  or  may  not  be  de¬ 
tectable,  depending  whether  Rx  <  RA  for  which  it  can  be  detected  (see 
above),  or  whether  R1  >  RA  for  which  it  cannot.  If  the  detector  does  not 
have  sufficient  sensitivity  to  detect  the  tank  at  range  Rlt  then  it  remains 
undetected  at  shorter  ranges,  until  it  emerges  from  the  fog  all  at  once  at  a 
relatively  large  image  size,  again  in  accord  with  common  experience.  When 
target  recognition  is  dominated  by  atmospheric  attenuation  rather  than  pix¬ 
els  on  target,  the  recognition  range  increases  with  detector  sensitivity  as 
=  -RaMC/NEP),  where  C  is  a  constant  determined  by  the  target  emis¬ 
sion  strength  and  camera  optics.  This  expression  can  be  more  complex  if 
forward  scattering  is  important.  The  recognition  range  Rr  is  determined 
primarily  by  atmospheric  absorption  through  RA ,  and  secondarily  by  detec¬ 
tor  sensitivity  through  the  logarithmic  dependence  on  NEP. 

Connections  between  target  recognition  range  and  detector  sensitivity 
have  important  implications  for  the  choice  of  detector  technology.  High  sen¬ 
sitivity  is  desirable  for  long  recognition  ranges  under  poor  atmospheric  con¬ 
ditions,  and  for  cameras  with  small  numerical  apertures:  it  appears  that 
HgCdTe  will  remain  the  technology  of  choice  for  infrared  telescopes,  long 
range  sights,  and  other  demanding  applications.  Uncooled  arrays  have  be¬ 
come  competitive  with  HgCdTe  when  recognition  at  very  long  ranges  is  not 
required,  and  uncooled  technology  may  well  come  to  dominate  the  majority 
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of  infrared  camera  applications.  When  the  sensitivity  of  uncooled  arrays  is 
adequate  to  produce  a  good  image  with  available  optics,  the  lower  NEP  of 
HgCdTe  detectors  has  little  advantage  in  recognition  range:  Rr  is  nearly  in¬ 
dependent  of  NEP  in  a  clear  atmosphere,  and  increases  only  logarithmically 
with  NEP  with  atmospheric  attenuation.  High  spatial  uniformity  and  the 
lack  of  defects  can  be  more  important  to  target  recognition  than  NEP  under 
these  conditions.  It  is  here  that  QWIP  detectors  have  an  advantage,  despite 
their  low  quantum  efficiency.  QWIP  technology  is  particularly  well  suited  to 
complex  pixel  designs  for  multi-color  and  polarization  sensing  arrays. 


3  UNCOOLED  MICROBOLOMETER  AR¬ 
RAYS 


The  development  of  uncooled  microbolometer  arrays  started  as  a  se¬ 
cret  program,  which  was  later  declassified  and  transferred  to  commercial 
production.  The  idea  is  simple:  to  use  modern  lithography  and  processing 
techniques  to  fabricate  large  arrays  of  small  temperature  sensing  bolometers. 
Two  approaches  were  originally  developed:  bolometer  arrays  with  resistive 
response  (Honeywell)  and  arrays  with  ferroelectric  response  (Texas  Instru¬ 
ments).  Figure  2  is  a  drawing  of  one  pixel  of  the  original  Honeywell  concept 
-  pixels  of  similar  geometry  are  now  made  in  a  monolithic  process  by  Texas 
Instruments  using  ferroelectric  sensors.  As  shown,  the  temperature  sensing 
element  is  thermally  isolated  above  the  readout  circuitry  by  using  long  thin 
legs,  which  also  provide  electrical  connections.  This  arrangement  provides  a 
high  fill  factor.  The  entire  array  is  made  in  a  monolithic  process,  based  on 
standard  Si  processing  with  the  addition  of  materials  needed  for  resistive  or 
ferroelectric  sensing  of  the  temperature. 


Circuitry 

Figure  2:  Diagram  of  a  bolometer  pixel  from  an  uncooled  bolometer  array. 
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Of  all  the  types  of  infrared  detectors,  uncooled  microbolometer  arrays 
are  likely  to  have  the  largest  impact  on  Army  operations.  That  is  because 
microbolometer  arrays  combine  very  good  quality  imagery  with  the  low  cost 
and  ease  of  use  associated  with  uncooled  detectors.  Although  their  sensi¬ 
tivity  does  not  approach  that  of  the  best  HgCdTe  detector  arrays,  uncooled 
bolometer  arrays  provide  a  level  of  performance  that  is  adequate  for  perhaps  a 
majority  of  applications.  Because  microbolometer  arrays  are  available  inter¬ 
nationally  on  a  commercial  basis,  one  can  no  longer  assume  that  an  adversary 
lacks  good  night  vision  capability. 


Although  uncooled  bolometers  axe  commonly  considered  to  be  insensi¬ 
tive  compared  with  cooled  quantum  detectors  such  as  HgCdTe,  they  can  pro¬ 
vide  surprisingly  good  performance.  As  detailed  in  Appendix  B.l,  the  square 
of  the  infrared  noise  equivalent  power  (NEP)  for  a  bolometer  is  (6P?h)  = 
4kBT2GAf,  where  kB  is  Boltzmann’s  constant,  G  is  the  thermal  conductance 
between  the  bolometer  element  and  the  substrate,  and  A /  is  the  bandwidth. 
For  quantum  detectors  the  square  of  the  NEP  is  (8P2)  =  2P1huAf  with  P7 
the  absorbed  infrared  power,  and  hv  the  photon  energy.  The  ratio  is 


(Of?)  _  1  /AT,-,  (  hv  \ 
(SP&)  2  \  T  )  \kBT) 


(3-1) 


where  ATy  is  the  temperature  increase  in  the  bolometer  produced  by  illumi¬ 
nation.  Comparing  the  two  types  of  detector,  we  find  that  the  square  of  the 
NEP  is  larger  for  quantum  detectors  by  the  ratio  of  the  photon  energy  to  the 
thermal  energy,  hu/kBT  ~  4  for  LWIR  detectors,  but  smaller  by  the  factor 
AT7/T.  Better  bolometer  performance  is  obtained  with  larger  A P7,  which 
can  be  increased  with  faster  optics  and  lower  thermal  conductance.  For  cur¬ 
rent  microbolometer  arrays  ATyl  K,  giving  SPy/5Pth  ~  0.1  and  an  ideal 
limiting  performance  NETD  ~  10  mK  for  uncooled  bolometers  vs.  NETD  ~ 
1  mK  for  quantum  detectors.  The  actual  performance  achieved  is  somewhat 
less  than  NETD  ~  75  mK,  due  to  excess  noise  in  the  bolometer  element  and 
readout  circuitry. 
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Areas  for  improvement  in  microbolometer  performance  are  higher  bolome¬ 
ter  responsitivity  (for  lower  noise),  better  thermal  isolation  (for  higher  in¬ 
duced  temperature  rise),  and  smaller  heat  capacity  (for  faster  response). 

High  responsitivity  materials  are  needed  for  bolometers  to  provide  a 
large  signal  and  large  signal-to-noise  ratio.  For  resistive  bolometers  the  re¬ 
sponsitivity  is  characterized  by  the  dimensionless  quantity  a  =  ( T/R)dR/dT . 
The  ratio  of  intrinsic  thermal  fluctuation  noise  to  Johnson  noise  in  the  tem¬ 
perature  sensor  is 


where  IB  is  the  bias  current,  5R  is  the  resistance  change  due  to  thermal  fluc¬ 
tuations,  5Vj  is  the  Johnson  noise  voltage,  and  ATg  is  the  temperature  rise 
produced  by  the  bias  current.  From  this  expression  we  see  that  high  responsi¬ 
tivity  and  large  bias  currents  are  desirable  to  overcome  the  intrinsic  Johnson 
noise  of  resistive  sensor  elements.  Comparing  different  classes  of  resistive 
materials  (see  Appendix  B.l)  we  find  a  ~  1  for  metals,  a  rsj  —  Eact/kBT  for 
semiconductors  with  Eact  the  activation  energy,  and  a  —  (1/4)  (To/T)1/4 
for  variable  range  hopping  materials  like  VOx.  The  parameter  To  character¬ 
izes  variable  hopping  in  the  material  used.  Cooled  semiconductors  provide 
excellent  responsitivity,  but  the  best  response  for  uncooled  resistive  detectors 
is  obtained  with  materials  like  VOx  for  which  a  ~  6  .  Using  this  value  of 
responsitivity  we  find  that  significant  heating  by  the  bias  current  A TB  ~  8 K 
is  needed  to  avoid  domination  by  Johnson  noise.  To  avoid  excessive  power 
consumption  by  the  array,  the  bias  current  can  be  pulsed  during  readout. 

Materials  that  exhibit  a  phase  transformation  near  the  temperature  of 
operation  are  natural  candidates  to  provide  higher  thermal  responsitivity. 
The  approach  used  by  Texas  Instruments  relies  on  a  ferroelectric  phase  tran¬ 
sition  near  room  temperature  to  obtain  high  responsitivity.  Phase  transi¬ 
tions  in  resistive  sensors  are  good  candidates  to  improve  performance.  The 
superconducting  transition  in  high  Tc  materials  results  in  a  strong  thermal 
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response,  but  current  transition  temperatures  require  cooling.  New  “colos¬ 
sal”  magnetoresistive  materials  developed  for  magnetoresistive  heads  have  an 
unusually  large  thermal  response  near  room  temperature,  as  shown  in  Fig¬ 
ure  3  (Venkatesan  et  al.,  University  of  Maryland).  This  material  (LCMO  - 
Lanthanam  Calcium  Manganese  Oxide)  shows  a  dimensionless  temperature 
coefficient  of  resistance  a  ~  40  at  T  -  275  K.  At  present  the  1/f  noise  level 
in  this  material  is  too  high  to  be  useful  in  detectors,  but  reductions  in  noise 
through  better  processing  may  lead  to  improved  sensors. 


Epitaxial  LCMO  FQm  :  Bolometric  Optical  Response 


Figure  3:  Resistive  transition  in  colossal  magnetoresistance  material  LCMO 
(Lanthanum  Calcium  Manganese  Oxide)  (T.  Venkatesan,  R.  Ramesh  and  M. 
Rajeswari,  Univerisity  of  Maryland). 


LWIR  infrared  cameras  based  on  uncooled  microbolometer  arrays  are 
now  commercially  available  from  a  number  of  companies  including  Santa 
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Barbara  Research  Center  (SBRC)  and  Texas  Instruments.  The  performance 
of  the  SBRC  camera  is  representative  of  the  better  units  -  it  is  based  on  a 
uncooled  resistive  320  x  240  pixel  microbolometer  array  with  response  in  the 
LWIR  band  (7  to  13.5  /im).  The  sensitivity  of  the  pixels  in  the  array  (NETD 
~  50  mK  with  f/1.0  optics)  and  the  uniformity  of  the  array  are  both  very 
good.  The  subjective  quality  of  video  images  taken  with  the  camera  is  high, 
better  than  first  generation  HgCdTe  imagers,  and  one  expects  that  uncooled 
LWIR  cameras  will  be  useful  for  a  wide  range  of  imaging  applications. 

Improvements  in  uncooled  microbolometer  arrays  are  underway  in  a 
number  of  areas.  Higher  responsitivity  is  possible  using  crystalline  VOx  in¬ 
stead  of  disordered  VOx  in  resistive  bolometers,  and  possibly  through  the  use 
of  colossal  magnetoresistive  materials.  Better  thermal  isolation  of  the  pix¬ 
els  can  be  achieved  through  narrower  and  longer  pixel  ‘legs’  leading  to  larger 
thermal  response  and  lower  noise.  The  ultimate  limit  is  determined  by  black- 
body  thermal  coupling  between  the  pixel  and  its  environment.  Larger  arrays 
(640x480  and  above)  with  smaller  pixels  (25x25  /im2)  are  under  develop¬ 
ment  to  increase  the  number  of  pixels  on  target.  The  combination  of  better 
responsitivity  and  better  thermal  isolation  should  give  improved  sensitivity 
NETD  ~  20  to  30  mK.  Improved  readout  IC’s  can  give  lower  noise  and 
better  uniformity.  Power  consumption  is  an  important  issue  for  uncooled 
cameras,  because  they  are  inherently  portable,  and  could  be  widely  used  in 
the  field.  The  focal  plane  array  consumes  ~300  mW,  while  the  rest  of  the 
electronics  and  display  take  5  W  to  10  W.  Lower  power  operation  should  be 
possible  through  the  use  of  low  power  electronics.  Low  cost  cameras  are  also 
under  development.  The  current  price  $10k  to  $20k  could  fall  to  ~$lk  if 
large  volume  orders  are  placed  for  mass  market  applications  such  as  those  in 
the  automotive  industry.  A  large  fraction  of  the  cost  and  weight  of  current 
uncooled  cameras  is  associated  with  the  Ge  lens.  Reflective  optics  based  on 
metal  coated  precision  molded  plastics  could  lead  to  large  cost  and  weight 
reductions. 
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In  summary,  uncooled  LWIR  cameras  are  becoming  increasingly  com¬ 
petitive  with  HgCdTe  cameras  for  a  wide  range  of  applications.  The  at¬ 
tractive  qualities  of  portability  and  low  cost  could  lead  to  uncooled  cam¬ 
eras  dominating  the  market  for  LWIR  cameras,  leaving  only  the  high  end  to 
HgCdTe  devices.  Once  uncooled  cameras  are  widely  available  internation¬ 
ally,  one  should  no  longer  assume  that  an  adversary  lacks  good  night  vision 
equipment. 
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4  QUANTUM  WELL  INFRARED 
PHOTODETECTORS 


Quantum  well  infrared  photodetectors  (QWIPs)  are  man-made  extrinsic 
photoconductors  in  which  quantum  wells  replace  impurity  atoms.  In  a  con¬ 
ventional  extrinsic  photoconductor,  free  carriers  excited  from  dopant  atoms 
by  infrared  light  provide  the  signal;  in  QWIPs,  free  carriers  excited  from 
doped  quantum  wells  provide  the  signal.  The  ability  to  vary  the  binding 
energy  of  electrons  in  QWIPs  to  match  the  desired  IR  response  by  chang¬ 
ing  quantum  well  depth  and  width  is  an  important  advantage.  Most  QWIP 
devices  are  made  from  GaAs/AlGaAs  heterostructures  grown  by  molecular 
beam  epitaxy  (MBE)  or  metal  organic  chemical  vapor  deposition  (MOCVD). 

Figure  4  is  a  schematic  diagram  of  a  multiple  quantum  well  QWIP 
device.  The  well  depth  and  the  widths  of  the  wells  and  barriers  are  adjusted 
so  that  the  ground  state  of  electrons  in  the  well  is  bound  and  the  first  excited 


CONDUCTION  BAND 


Figure  4:  Schematic  diagram  showing  the  principle  of  operation  of  quantum 
well  infrared  photodetector  (QWIP)  devices  with  inset  of  scanned  electron 
microscope  photograph  of  GaAs/AlGaAs  heterostructure. 
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state  lies  near  the  top  of  the  well  as  indicated.  The  excited  states  overlap 
significantly  between  wells  and  form  the  analog  of  a  conduction  band  in 
extrinsic  semiconductors.  Absorbed  infrared  photons  promote  electrons  from 
the  bound  state  to  the  conduction  band  where  they  are  swept  away  by  an 
applied  electric  field  to  produce  the  photocurrent. 

The  fundamental  advantage  of  QWIP  devices  is  the  ability  to  tailor 
their  infrared  response,  and  they  are  naturally  suited  to  the  fabrication  of 
multi-color  sensors.  QWIPs  based  on  GaAs/AlGaAs  heterostructures  also 
benefit  from  the  high  quality  and  uniformity  of  these  semiconductor  struc¬ 
tures,  which  are  made  using  a  mature  technology.  As  discussed  below,  QWIP 
devices  are  quite  sensitive  to  the  polarization  of  infrared  photons  and  are  nat¬ 
ural  candidates  for  polarization  sensitive  sensors. 

The  primary  disadvantages  of  QWIP  devices  are  low  quantum  efficiency 
and  low  operating  temperature.  Quantum  efficiency  is  limited  by  a  number  of 
factors.  Because  each  well  contains  few  electrons  compared  with  intrinsic  ma¬ 
terial  and  absorbs  infrared  radiation  only  weakly,  a  stack  of  many  (typically 
~  100)  quantum  wells  is  required  in  order  to  obtain  sufficient  absorption. 
In  addition,  the  photoconductive  gain  (see  Appendix  B.2)  of  QWIPs  is  lim¬ 
ited  to  values  less  than  unity  by  the  fact  that  free  carriers  are  recaptured 
before  travelling  the  full  width  of  the  stack.  Quantum  mechanical  selection 
rules  limit  infrared  absorption  to  radiation  with  electric  field  polarized  per¬ 
pendicular  to  the  quantum  well.  An  optical  coupling  structure  (see  below) 
is  required  to  rotate  the  polarization  of  incoming  radiation  into  the  correct 
orientation  for  absorption. 

At  sufficiently  low  temperatures  QWIPs  can  give  excellent  sensitivity. 
However,  the  required  operating  temperatures,  ~70  K  and  below,  are  low 
relative  to  those  for  intrinsic  photoconductors  and  bolometers.  As  described 
in  Appendix  B.2,  extrinsic  photoconductors  require  relatively  low  operating 
temperature  for  fundamental  reasons.  Thermal  ionization  of  carriers  bound 
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in  dopant  atoms  (or  quantum  wells)  occurs  at  temperatures  ksT  well  be¬ 
low  the  binding  energy  (and  infrared  photon  energy),  because  the  entropy 
of  free  carriers  is  much  greater.  This  effect  is  compounded  by  the  fact  that 
the  infrared  induced  charge  carrier  density  -  the  signal  -  is  limited  by  short 
free  carrier  lifetimes  in  QWIP  devices.  Free-to-bound  transitions  in  quantum 
wells  are  very  rapid  due  to  the  nesting  of  electron  subbands  in  momentum 
space  with  free  carrier  lifetimes  typically  r  ~  10  ps,  whereas  the  minor¬ 
ity  carrier  lifetime  in  HgCdTe  materials  is  much  longer,  typically  r  ~  lfxs. 
Miniature  coolers  which  reach  the  required  temperatures  are  available,  but 
the  added  stress  limits  their  lifetime  to  ~2000  hours. 

Optical  couplers  are  required  for  all  QWIP  devices,  because  quantum 
well  transitions  are  only  excited  by  electric  fields  perpendicular  to  the  plane  of 
the  well.  A  coupler  is  needed  to  bend  the  infrared  rays  and  rotate  the  polar¬ 
ization.  Two  examples  of  optical  couplers  are  shown  in  Figure  5:  a  corrugated 
QWIP  developed  by  Princeton  University  and  the  Army  Research  Labora¬ 
tory  (ARL)  and  a  pseudorandom  coupler  developed  at  the  Jet  Propulsion 
Laboratory  (JPL).  The  corrugated  QWIP  relies  on  total  internal  reflection: 


Figure  5:  Optical  couplers  for  QWIP  devices:  (left)  corrugated  QWIP  devel¬ 
oped  by  Princeton  University  and  Army  Research  Laboratory;  (right)  pseu¬ 
dorandom  coupler  developed  by  the  Jet  Propulsion  Laboratory. 


infrared  radiation  incident  from  the  bottom  is  trapped  by  the  corrugations 
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on  the  top  of  the  device,  and  its  polarization  is  rotated.  The  parallel  roof 
shaped  structure  shown  is  polarization  sensitive  with  an  extinction  ratio  ~4. 
In  order  to  avoid  polarization  sensitivity,  JPL  has  built  couplers  with  pseu¬ 
dorandom  corrugations  as  shown  in  Figure  5.  The  complex  pattern  rotates 
incident  rays  to  minimize  sensitivity  to  polarization. 

QWIP  cameras  have  been  developed  by  a  number  of  makers.  Although 
the  requirement  for  a  cooler  tends  to  increase  size,  weight  and  cost,  JPL  has 
developed  miniature  hand-held  QWIP  cameras  with  very  good  performance. 
Lockheed  Martin  has  demonstrated  a  video  camera  using  a  256x256  pixel 
QWIP  device  operating  at  60  K.  With  f/1.7  optics  and  a  2.27°  x  2.27°  field 
of  view  this  camera  achieves  a  NETD  varying  from  ~7  mK  to  35  mK  with 
spatial  frequency.  The  subjective  quality  of  the  images  is  excellent,  due  in 
large  part  to  the  high  spatial  uniformity  of  the  array. 

An  important  advantage  of  QWIP  devices  is  that  they  are  based  on  a 
mature  GaAs/AlGaAs  technology  which  produces  high  quality  material  and 
permits  the  growth  of  complex  layered  structures.  As  a  result,  the  uniformity 
of  QWIP  devices  is  very  high,  and  their  performance  approaches  theoreti¬ 
cal  limits.  Their  primary  disadvantages  are  low  quantum  efficiency,  and  the 
need  for  cooling  to  relatively  low  temperatures.  Consequently  their  sensi¬ 
tivity  is  limited.  Nonetheless  the  subjective  quality  of  QWIP  images  can 
be  much  better  than  that  for  other  infrared  cameras  with  higher  sensitivity 
but  lower  uniformity.  The  trade-off  between  sensitivity  and  uniformity  is 
an  interesting  issue  for  pattern  recognition,  as  described  above.  Despite  an 
uncertain  past,  polarization  of  infrared  radiation  may  prove  to  be  a  useful 
discriminant  for  recognition  of  metal  targets  in  organic  clutter.  QWIP  detec¬ 
tor  arrays  are  natural  candidates  for  this  application.  QWIP  detector  arrays 
are  also  naturally  suited  to  “multi-color”  infrared  detectors  arrays,  because 
their  wavelength  sensitivity  can  be  continuously  adjusted  during  growth. 
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5  HgCdTe  PHOTODETECTOR  ARRAYS 


HgCdTe  photodetectors  have  fundamental  advantages  for  infrared  sen¬ 
sors  (detailed  in  Appendix  B.2).  Because  they  are  intrinsic  photoconductors 
in  which  infrared  photons  excite  transitions  from  the  valence  to  the  con¬ 
duction  band,  they  have  high  quantum  efficiency  and  low  thermally  excited 
carrier  density  (relative  to  extrinsic  photoconductors).  Their  band  gap  can 
be  adjusted  by  varying  the  relative  proportion  of  Hg  and  Cd  to  produce  peak 
sensitivity  from  the  LWIR  through  the  MWIR  bands.  These  fundamental 
advantages  have  led  to  a  considerable  investment  in  HgCdTe  technology  over 
the  years.  However,  the  actual  performance  of  LWIR  HgCdTe  detectors  has 
not  reached  these  fundamental  limits  for  a  variety  of  quite  technical  reasons. 
For  alloys  appropriate  to  the  LWIR  band,  HgCdTe  is  not  highly  stable,  and 
the  quality  of  devices  is  very  sensitive  to  the  details  of  the  fabrication  process. 

An  important  development  between  this  JASON  study  and  the  previ¬ 
ous  one  is  the  production  of  large  LWIR  staring  arrays  made  from  HgCdTe. 
Both  Santa  Barbara  Research  Center  and  Texas  Instruments  have  produced 
arrays  of  size  256x256  pixels  and  640x480  pixels  which  give  excellent  im¬ 
age  quality.  Reliability  has  been  improved  by  providing  strain  relief  in  the 
connections  between  the  HgCdTe  array  and  the  Si  readout  integrated  circuit 
(ROIC).  SBRC  attacks  this  problem  by  backing  the  Si  readout  with  propri¬ 
etary  materials  to  match  the  thermal  expansion  of  the  HgCdTe  array.  Texas 
Instruments  epoxies  thinned  HgCdTe  arrays  to  Si  readout  ICs  so  that  the 
differential  strain  is  taken  up  by  the  HgCdTe  array. 

These  new  HgCdTe  LWIR  staring  arrays  provide  excellent  sensitivity 
along  with  high  spatial  resolution,  and  they  are  natural  choices  for  infrared 
telescopes  for  which  atmospheric  absorption,  and  scattering  and  absorption 
by  dust  and  aerosols  is  an  issue.  Once  infrared  capability  becomes  common, 
HgCdTe  arrays  will  provide  better  sensitivity  than  uncooled  bolometer  arrays 
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or  QWIPs.  However,  the  market  for  high  performance  HgCdTe  arrays  is 
small,  and  it  seems  that  their  price  will  remain  high,  especially  compared 
with  uncooled  detector  arrays.  There  is  a  significant  economic  issue  as  to 
who  will  pay  to  bring  large  HgCdTe  arrays  into  production. 
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6  DEVELOPMENTS 


This  section  describes  a  variety  of  developments  which  are  not  limited 
to  any  one  detector  type:  ‘two  color’  infrared  imagers,  improved  readout 
integrated  circuits,  improved  strain  relief,  and  improved  materials  for  ther¬ 
moelectric  coolers. 

The  use  of  ‘color’  in  infrared  images  -  the  presentation  of  multiple  wave¬ 
length  bands  in  one  image  -  promises  to  improve  discrimination  of  objects 
just  as  it  does  at  visible  wavelengths.  Different  materials  have  different  emis- 
sivity,  and  many  materials  have  characteristic  absorption  features  in  the 
LWIR  and  MWIR  bands.  Comparison  of  LWIR  and  MWIR  emission  can 
also  improve  the  accuracy  of  inferred  temperatures.  The  complexity  of  mul¬ 
tiple  color  infrared  imagers  using  separate  detector  arrays  for  different  ‘colors’ 
has  limited  their  application  in  the  past. 

Infrared  detector  arrays  in  which  two  colors  are  integrated  in  a  single 
pixel  are  under  development  for  QWIP  and  HgCdTe  technologies.  Both  are 
well  suited  to  multiple  color  detectors,  because  their  peak  sensitivity  can  be 
adjusted.  In  principle,  two  color  bolometer  arrays  could  also  be  fabricated 
by  using  wavelength  selective  coatings.  Figure  6  illustrates  two  approaches 
to  the  fabrication  of  two-color  pixels. 

In  a  QWIP  device  the  wavelength  sensitivity  can  be  changed  by  adjust¬ 
ing  the  quantum  well  depth  and/or  width  during  growth  of  the  GaAs/AlGaAs 
heterostructure.  Figure  6  (top)  illustrates  a  two  color  detector  pixel  in  which 
two  QWIP  photoconductors  are  grown  in  series.  Selective  patterning  and 
etching  techniques  are  used  to  contact  the  two  detectors  separately  as  indi¬ 
cated.  For  a  HgCdTe  device  the  wavelength  sensitivity  is  adjusted  by  chang¬ 
ing  the  alloy  concentrations.  Figure  6  (bottom)  shows  a  hetero junction  pixel 
that  provides  simultaneous  MWIR  and  LWIR  sensitivity. 
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Figure  6:  Two  color  infrared  detector  pixels:  (top)  QWIP  two  color  pixel 
consisting  of  stacked  MWIR  and  LWIR  multiple  quantum  wells  grown  in 
a  single  GaAs/AlGaAs  heterostructure;  (bottom)  HgCdTe  heterojunction 
detector  providing  simultaneous  MWIR  and  LWIR  band  two  color  capability 
with  spectral  response  shown. 
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Readout  integrated  circuits  (ROICs)  have  improved  significantly  over 
the  past  five  years.  Silicon  CMOS  circuitry  has  generally  replaced  silicon 
charged  coupled  devices  (CCDs)  enabling  the  infrared  industry  to  use  stan¬ 
dard  foundry  facilities.  The  use  of  CMOS  lowers  cost  and  enables  the  design 
of  higher  capability  readouts.  Current  CMOS  designs  typically  place  the 
amplifier  front  end  and  a  switch  under  each  pixel.  A  row  of  preamps  and 
analog-to-digital  converters  (ADCs)  is  placed  along  one  side  of  the  chip,  and 
the  array  is  read  out  by  sweeping  the  active  pixels.  Unless  suitable  integra¬ 
tion  can  be  placed  under  each  pixel,  this  multiplexing  technique  reduces  the 
signal-to-noise  ratio  due  to  inefficient  time  averaging.  Economic  factors  often 
limit  infrared  detector  array  developers  to  the  use  of  old  CMOS  processes 
and  relatively  small  wafer  sizes. 

Considerable  improvement  would  be  possible  if  current  CMOS  process¬ 
ing  were  used  for  infrared  readout  ICs.  For  example,  with  0.25  /im  CMOS 
one  could  place  an  entire  16-bit  ADC  under  each  pixel.  With  such  an  ar¬ 
rangement  each  pixel  could  be  continuously  read  out  for  improved  signal  to 
noise.  One  could  also  place  considerable  processing  power  on  the  readout 
IC  to  use  for  signal  processing  and  feature  recognition.  The  availability  of 
these  new  capabilities  could  have  substantial  implications  for  new  readout 
IC  designs. 

The  size  of  current  readout  ICs  limit  the  size  of  infrared  arrays.  For 
example,  a  1000x1000  pixel  array  with  25x25  fxm2  pixels  has  overall  dimen¬ 
sions  of  2.5 x 2.5  cm2  or  lxl  in2.  This  is  about  the  limit  of  current  ROIC 
technology,  but  could  be  improved  substantially  by  using  modern  processes, 
to  produce  arrays  with  more  or  larger  pixels. 

Strain  relief  is  important  for  the  reliability  of  cooled  infrared  detector 
arrays.  Differential  thermal  contraction  between  the  detector  array  and  the 
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Si  readout  IC  places  stress  on  the  interconnects  which  can  lead  to  premature 
failure.  A  variety  of  approaches  have  been  developed  over  the  past  five  years 
to  address  this  issue,  as  illustrated  in  Figure  7. 


Interconnect 


Figure  7:  Approaches  to  thermal  strain  relief:  (left)  thinned  QWIP  arrays 
indium  bump  bonded  to  Si  readout  IC;  (right)  thinned  HgCdTe  array  epoxy 
bonded  to  Si  readout  IC. 
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A  FUNDAMENTALS  OF  INFRARED  DE¬ 
TECTOR  PERFORMANCE 


In  this  section  we  briefly  discuss  fundamental  considerations  which  limit 
infrared  detector  performance.  Subsection  A.l  on  blackbody  radiation  presents 
the  spectrum  and  noise  limits;  and  Subsection  A.2  discusses  the  physical  con¬ 
siderations  which  determine  the  detector  operating  temperature.  In  both  we 
emphasize  simple  physical  models  which  can  be  used  to  compare  different 
detector  types. 

A.l  Blackbody  Radiation 


Spectrum 


For  detectors  of  thermal  infrared  radiation,  blackbody  radiation  is  both 
the  source  of  the  signal  and  an  important  source  of  noise.  The  spectrum 
of  blackbody  radiation  is  easily  derived  from  the  statistics  of  photons  in 
equilibrium  with  a  blackbody  of  temperature  T.  The  volume  density  of 
photons  dN-y/d\  per  wavelength  interval  dX  is  simply  the  product  of  the 
density  of  photon  states  87r/A4  and  the  average  occupation  of  each  state 
n(hc/ X)  =  l/[exp(hc/XksT)  —  1],  the  Bose-Einstein  distribution  function: 


dNy  /87t\  1 

~dX  =  v  AV  ehc/xkzT  -  1  ’ 


(A-l) 


where  he/ X  is  the  photon  energy,  with  h  Planck’s  constant  and  c  the  velocity 
of  light,  and  ks  is  Boltzmann’s  constant. 


The  blackbody  photon  flux  dl^/dX  incident  on  one  side  of  a  flat  detec¬ 
tor  of  area  A  collected  from  a  solid  angle  f2  symmetric  about  the  direction 
perpendicular  to  the  detector  is: 


d/7 

dX 


=  A 


C\  dNry 

V  ~dX’ 


(A-2) 
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which  reduces  to: 

d/7  .  r /0\  /fl\2l  c  1  ,  x 

d\  V  7T  /  V27T/  A4  ehc/\kBT  _  J '  (A3) 

Plots  of  blackbody  photon  flux  spectra  are  shown  in  Figure  A-l  for  a 
series  of  temperatures  T  =  240K  to  340K  assuming  a  solid  angle  Q  =  2tt. 
As  shown,  the  peak  photon  flux  occurs  near  A  =  10/um,  and  falls  off  rapidly 
at  shorter  wavelengths.  In  order  to  detect  a  blackbody  signal  from  a  warm 
object,  one  would  like  to  absorb  as  much  of  the  radiation  as  possible,  and 
response  in  the  LWIR  atmospheric  transmission  band  near  lOfim  is  desirable. 
As  shown,  the  photon  flux  in  the  5//m  MWIR  band  is  appreciable  for  tem¬ 
peratures  T  >  300K,  but  falls  exponentially  fast  at  lower  temperatures.  This 
lack  of  response  of  MWIR  detectors  below  T  ~  300K  can  cause  difficulties 
in  imaging  relatively  cool  objects  and  backgrounds. 


Blackbody  Photon  Flux 


Wavelength  (micron) 

Figure  A-l:  Blackbody  photon  flux  vs.  wavelength  at  various  temperatures 
indicated  (see  text) 
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Fluctuations  in  Photon  Flux 


Photons  are  continually  emitted  and  absorbed  by  a  blackbody,  leading 
to  thermal  fluctuations  5n  in  the  occupation  n(hc/X)  of  each  quantum  state 
of  the  photon  gas  according  to  Boson  statistics: 

(Sn2)  =  n(l  +  n).  (A-4) 

Near  the  peak  of  the  blackbody  spectrum,  the  occupation  n  <  1,  the  photon 
gas  is  nondegenerate,  and  the  fluctuations  ( 5n 2)  =  n  obey  simple  counting 
statistics.  At  longer  wavelengths  (A  »  50/im  for  T  =  300K)  the  photon  gas 
becomes  degenerate  (n  >  1)  and  fluctuations  in  occupation  are  much  larger 
(Sn2)  =  n2;  this  regime  corresponds  to  fluctuating  speckle  patterns  with  a 
high  degree  of  coherence. 

The  total  number  of  photons  which  arrive  at  a  detector  over  an  integra¬ 
tion  time  At  is: 

iV7  =  ^~AXAt,  (A-5) 

(LA 

where  AA  is  the  wavelength  interval  detected.  Assuming  that  the  photon 
gas  is  nondegenerate  at  the  wavelengths  of  interest,  the  fluctuations  SN7  in 
number  of  photons  are  given  by: 

(SN2)  =  Ny.  (A-6) 

These  fluctuations  correspond  to  shot  noise  in  the  arrival  rate  J7  of  pho¬ 
tons.  For  a  given  bandwidth  Af  in  Hz,  fluctuations  in  the  photon  current  at 
frequencies  much  lower  than  the  arrival  rate  are  given  by: 

<«?>  -  21-Af-  (A-7) 

Generation-Recombination  Noise 

Now  we  consider  the  effect  of  the  incident  photon  flux  on  an  infrared 
detector,  which  we  take  to  be  a  photoconductor.  Each  incident  photon  cre¬ 
ates  a  photoexcited  carrier  with  probability  r],  the  quantum  efficiency.  For 


29 


simplicity  we  assume  that  only  one  type  of  excited  carrier  (electrons  or  holes) 
is  mobile.  After  creation,  the  photoexcited  carrier  drifts  in  the  electric  field 
applied  between  the  contacts  of  the  detector  for  a  distance  before  it  recom¬ 
bines.  In  doing  so,  it  transfers  a  charge  q  =  ( T0/rt)e  between  the  contacts 
and  through  the  external  circuit,  where  ra  is  the  carrier  lifetime,  rt  is  the 
transit  time  across  the  full  width  of  the  detector,  and  the  ratio  G  =  (r0/rf) 
is  the  photoconductive  gain.  Note  that  the  average  transferred  charge  q  is 
a  macroscopic  quantity  which  can  have  a  continuous  range  of  values  from 
q  <  e  to  q  >  e  determined  by  the  photoconductive  gain  (see  below).  The 
total  number  of  carriers  transferred  to  the  external  circuit  in  a  time  interval 
At  is: 

N  =  r/GNy,  (A-8) 

with  corresponding  photocurrent: 

I  =  rjGIy.  (A-9) 


Blackbody  radiation  incident  on  a  photoconductor  creates  carriers  ac¬ 
cording  to  the  process  described  above.  Fluctuations  in  the  number  of  charges 
transferred  to  the  circuit  in  a  time  At  are  determined  by  the  statistics  of  car¬ 
rier  generation  and  recombination,  as  well  as  fluctuations  in  the  photon  flux: 


5N2 

\  AT2 


>  -  <-^r>  +  <g5->  +  <-jy5">- 


(A-10) 


Here  ( Sr/ 2)  and  (5G2)  are  fluctuations  in  the  quantum  efficiency  77  and  pho¬ 
toconductive  gain  G  determined  by  averaging  over  the  arrival  of  Ny  photons. 
It  is  a  straightforward  exercise  in  statistics  to  show  that: 


-  T’  “d  (5G2)  =  S'  (A-n) 

Combining  Equations  (A-6),  (A-10),  and  (A-ll)  we  find: 

,6N\  (1  -  77)  ,  1  1  2 

(N-S  vNy  +  VNy  +  Ny~  VNy  (A'12) 


30 


Thus  fluctuations  ( 8N2)/N 2  in  the  number  of  carriers  transferred  to  the 
external  circuit  are  determined  by  the  number  rjNy  of  photoexcited  carriers, 
and  are  independent  of  the  value  G  of  the  photoconductive  gain.  The  factor 
of  2  in  Equation  (A-13)  arises  because  both  the  generation  and  recombination 
of  photocarriers  produce  statistical  noise.  The  photoconductive  gain  G  does 
not  introduce  additional  noise,  but  acts  in  the  same  way  as  amplifier  gain. 
The  corresponding  generation-recombination  noise  current  in  the  external 
circuit  at  frequencies  /  <C  1/27 rr0  is: 

(SI2)  =  AelGAf.  (A-13) 

We  assume  for  the  moment  that  the  thermal  generation  of  carriers  is  negligi¬ 
ble,  so  that  I  in  Equation  (A-13)  is  the  photocurrent  produced  by  blackbody 
background  radiation.  In  photodiodes,  no  statistical  noise  is  produced  by  re¬ 
combination,  because  excited  photocarriers  are  swept  out  of  the  active  region 
of  the  device  before  they  recombine;  for  this  case  G  =  1  and  the  noise  (SI2)  is 
a  factor  of  two  smaller  than  for  Equation  (A-13).  As  for  quantum  detectors, 
bolometers  are  also  ultimately  limited  by  photon  shot  noise,  although  this 
limit  is  not  commonly  approached  in  practice. 


The  fundamental  limit  to  the  performance  of  an  ideal  infrared  detector 
is  given  by  the  generation-recombination  noise  associated  with  background 
blackbody  radiation.  This  limit  is  usually  expressed  in  terms  of  the  noise 
equivalent  photon  power  (NEP)  incident  on  the  detector: 


Inserting  the  blackbody  photocurrent  fluctuations  (SI2)  due  to  generation- 
recombination  noise  in  photoconductors  from  Equation  (A-13)  we  find: 


1 

g/ic/AfcsT _  j 


1/2 

(A-15) 


For  photodiodes  the  noise  equivalent  power  due  to  photon  shot  noise  is 
smaller  than  Equation  (A-15)  by  a  factor  l/y/2.  These  expressions  are  valid 
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only  for  hc/\kBT  >  1,  because  we  have  assumed  that  the  photon  gas  is  non¬ 
degenerate  (incoherent).  This  assumption  is  justified  near  the  peak  of  the 
blackbody  spectrum  and  at  shorter  wavelengths.  For  blackbody  sources  the 
signal-to-noise  ratio  is: 


C/AT  —  (^C/A)J7 

/  NEP  ' 


(A-16) 


where  /7  is  the  blackbody  photon  current  from  Equation  (A-3).  Combining 
expressions  we  find  for  photoconductors: 
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(A-17) 


For  photodiodes  the  signal-to-noise  ratio  due  to  photon  shot  noise  is  larger 
than  Equation(A-17)  by  a  factor  \f2.  The  signal-to-noise  ratio  is  maximized 
for  thermal  blackbody  sources  by  maximizing  the  number  of  detected  black¬ 
body  photons  through  increases  in  detector  area  A,  solid  angle  ft,  spectral 
width  AA,  and  the  quantum  efficiency  77,  and  operation  at  wavelengths  A 
near  the  peak  of  the  blackbody  photon  flux  spectrum,  shown  in  Figure  A-l. 
Except  for  increasing  the  quantum  efficiency  77,  these  trends  also  increase  the 
NEP  and  would  be  undesirable  for  non-thermal  sources. 


Another  standard  measure  of  infrared  noise  performance  is  the  specific 
detectivity: 


D * 


/Al/2Afl/2\ 

\  NEP  ) ' 


(A- 18) 


The  background-limited  infrared  performance  (BLIP)  limit  from  Equation 
(A-15)  for  photoconductors  is: 


D* 


©-(DT1^-1))172-  (A-19) 


For  photodiodes  the  specific  detectivity  D*  due  to  photon  shot  noise  is  larger 
than  Equation(A-19)  by  a  factor  y/2.  For  T  =  300K,  ft  =  2ir,  77  =  1,  A  = 
10/un,  and  AA  =  1/zm,  the  BLIP  detectivity  for  photoconductors  from  Equa¬ 
tion  (A- 19)  is 

D*  =  6.43  x  1010cmHz1/2/IF.  (A-20) 
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Note  that  the  use  of  the  specific  detectivity  D*  to  characterize  detectors  of 
thermal  sources  can  be  misleading,  because  lower  NEP  is  not  always  equiva¬ 
lent  to  higher  signal  to  noise  ratio,  as  discussed  above.  For  example,  detectors 
which  have  smaller  solid  angles  Cl  or  do  not  respond  at  wavelengths  near  the 
peak  in  the  blackbody  spectrum  can  have  higher  detectivity  than  Equation 
(A-20),  but  have  lower  signal-to- noise  ratios  for  thermal  sources,  because 
they  detect  fewer  photons. 

A  measure  of  noise  performance  of  detector  systems  is  the  noise-equivalent 
temperature  difference  (NETD),  approximately  the  smallest  temperature 
change  which  can  be  detected  in  a  source  of  thermal  radiation.  The  NETD 
is  particularly  useful,  because  it  includes  the  influence  of  all  relevant  param¬ 
eters  and  is  an  unambiguous  measure  of  the  performance  of  a  given  system. 
The  NETD  is  defined  as  the  temperature  difference  of  a  blackbody  source 
at  a  specified  temperature  (e.g.  T  =  300  K)  which  corresponds  to  the  same 
change  in  detected  signal  as  the  system  noise. 

A. 2  Temperature  of  Operation 

Most  infrared  detectors  with  sensitivity  in  the  LWIR  band  require  cool¬ 
ing  in  order  to  reduce  the  number  of  thermally  excited  carriers.  The  tem¬ 
perature  of  operation  required  for  a  certain  level  of  performance  is  a  critical 
parameter,  because  the  cooling  system  often  dominates  the  size,  weight  and 
reliability  of  the  detector  system.  Temperatures  above  the  77K  boiling  point 
of  liquid  nitrogen  are  desirable,  because  they  can  be  reached  using  nitrogen- 
based  coolers  which  are  standard  in  many  applications.  Lower  temperatures 
can  be  reached  using  other  gases,  at  the  expense  of  greater  weight  and  non¬ 
standard  operation. 

The  specific  factors  which  determine  the  temperature  of  operation  for 
a  given  detector  system  are  complex,  and  it  is  not  always  easy  to  compare 
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the  promise  of  different  approaches.  In  this  section  we  discuss  the  physical 
processes  which  determine  operation  temperature  from  a  general  point  of 
view,  and  present  a  number  of  simple  models.  Our  goal  is  to  obtain  an 
approximate  picture  which  is  valid  for  a  range  of  detector  designs,  rather 
than  a  detailed  description  of  any  given  detector. 

Quantum  Detectors 

Generation-recombination  (G-R)  noise  associated  with  the  production 
of  thermally-excited  carriers  sets  the  fundamental  limit  to  quantum  detector 
performance,  as  discussed  above  in  Section  2.  As  the  detector  is  cooled, 
the  G-R  noise  level  decreases  as  the  square  root  of  the  thermally-generated 
carrier  density  nth,  and  the  signal-to-noise  ratio  increases  proportionately. 
The  physical  quantity  which  determines  the  signal-to-noise  ratio  is  the  ratio 
ns/nth  of  the  carrier  concentration  ns  produced  by  the  infrared  signal  to  the 
thermally-generated  concentration  nth ■  A  primary  goal  in  detector  design  is 
to  increase  this  ratio  as  much  as  possible. 

In  this  subsection  we  consider  the  statistics  of  carrier  generation  and  re¬ 
combination  in  order  to  develop  principles  for  the  design  of  infrared  detectors 
with  high  operation  temperatures.  Using  simple  models,  we  show  that  the 
operating  temperature  is  maximized  by  maximizing  the  radiative  efficiency 
of  the  carriers  generated  by  infrared  radiation,  by  matching  the  density  of 
levels  between  which  radiative  transitions  take  place.  It  is  well  known  that 
matching  the  absorption  spectrum  of  the  detector  to  the  source  of  radiation 
also  maximizes  operating  temperature  by  reducing  the  thermal  generation 
rate  associated  with  states  which  absorb  blackbody  radiation  but  are  not 
excited  by  the  signal. 

For  simplicity  we  consider  a  two-level  model  of  a  quantum  detector. 
The  lower  level  represents  the  ground  state  of  electrons  in  the  system,  and  is 
assumed  to  be  heavily  populated  at  low  temperatures;  specific  examples  are 
a  donor  level  for  an  extrinsic  photoconductor,  or  a  bound  state  in  a  quantum 
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well.  The  upper  level  represents  the  excited  states  of  the  system  which  are  as¬ 
sumed  to  be  unpopulated  at  sufficiently  low  temperatures;  specific  examples 
are  the  conduction  bands  of  extrinsic  photoconductors,  or  of  quantum  well 
detectors.  We  assume  that  only  excited  carriers  are  mobile  and  contribute  to 
the  detected  signal;  the  free  carrier  concentration  is  labelled  n  and  we  refer  to 
these  as  electrons.  The  volume  density  of  quantum  states  in  the  lower  level 
is  labelled  Ni  and  that  in  the  upper  level  A^.  We  assume  that  a  fraction  of 
the  lower  states  with  density  No  are  empty  even  at  very  low  temperatures 
due  to  the  presence  of  occupied  traps  which  remove  carriers  from  the  system. 
For  example,  an  n-type  extrinsic  photoconductor  with  donor  and  acceptor 
concentrations  No  and  Na  would  have  No  =  Na,  Ni  =  No,  and  N2  =  Nc, 
where  Nc  is  the  effective  density  of  states  in  the  conduction  band,  given 
below. 

Radiative  Efficiency 

For  the  ideal  case  in  which  all  transitions  between  levels  1  and  2  occur 
via  radiative  processes,  the  rate  of  change  of  the  electron  density  n  is  given 
by: 

di~i 

—  =  ( JtVar  +  JsVar)(Ni  -  N0  -  n)  -  nvacr(N0  +  n).  (A-21) 

For  this  case  thermal  equilibrium  is  established  via  the  absorption  and  emis¬ 
sion  of  blackbody  radiation  at  the  detector  temperature  T,  and  no  other 
thermalization  processes  are  active.  The  first  thermal  generation  term  is  the 
rate  via  the  absorption  of  blackbody  radiation,  for  which  Jr  is  the  blackbody 
radiation  flux  and  oar  is  the  radiative  absorption  cross  section.  The  second 
generation  term  in  Equation(A-21)  is  the  rate  via  absorption  of  the  infrared 
signal  with  photon  flux  Js;  the  absorption  cross  section  for  the  signal  is  nor¬ 
mally  the  same  as  for  blackbody  radiation.  The  third  term  in  Equation(A-21) 
describes  the  radiative  recombination  of  excited  carriers  with  thermal  veloc¬ 
ity  v  and  capture  cross  section  a  a-.  If  the  density  of  excited  carriers  n  is 
much  smaller  than  the  trap  density  N0,  the  rate  equation  is  approximately 
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linear: 


an  77, 

—  =  (Jt<Jar  +  JsVarjiNx  -  N0)  -  —  (A-22) 

where  tt  =  l/(vcrcrNQ)  is  the  radiative  lifetime  of  excited  carriers.  The  solu¬ 
tions  to  Equation  (A-22)  are  simply  decaying  exponentials  with  characteristic 
time  tt. 


In  thermal  equilibrium  with  no  signal  present,  the  thermal  rate  of  gen¬ 
eration  of  carriers  just  balances  the  rate  of  recombination.  This  fact  can  be 
used  to  relate  the  cross  sections  for  radiative  absorption  and  capture.  In 
thermal  equilibrium  the  carrier  concentration  from  Equation  (A-22)  is: 


=  Jt<Jar(Ni  -  N0)tt. 


(A-23) 


The  thermal  carrier  concentration  nt  can  also  be  found  from  the  statistics  of 
electrons.  At  low  temperatures  for  which  relatively  few  carriers  are  excited 
(nt  <C  No)  we  have  approximately: 


nt  =  -^2(^1  Np  )_e~AE/kBT 

N0 

Combining  Equations  (A-23)  and  (A-24)  we  find  that  the 
sections  for  absorption  and  capture  are  proportional: 


(A-24) 
radiative  cross 


-  ( £)  ®  e_AE/t'T  - (x) (A-25> 

The  apparent  exponential  temperature  dependence  of  aar  in  Equation  (A- 
25)  is  cancelled  by  that  of  the  blackbody  flux  Jt,  and  the  cross  section  is 
independent  of  the  trap  density  Nq  as  is  required 


For  the  ideal  case  in  which  all  transitions  are  radiative,  the  ratio  ns/nt 
of  the  signal  to  the  thermally-generated  carrier  concentration  from  Equation 
(A-22)  is  simply: 


rh  _  Js 
nt  Jt' 


(A-26) 


Thus  the  flux  of  infrared  photons  from  the  signal  must  be  sufficiently  large 
compared  with  that  from  a  blackbody  at  the  detector  temperature  in  order 
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to  obtain  good  noise  performance.  However,  this  ideal  case  is  not  always 
obtained  in  practice,  because  thermalization  of  the  carrier  distribution  often 
occurs  via  non-radiative  mechanisms  such  as  phonon  absorption  and  emis¬ 
sion.  For  this  case  we  need  to  add  the  nonradiative  rates  of  carrier  generation 
and  recombination  to  the  rate  equation  for  n: 

%  =  {JtOar  +  JsVar  +  -N0)-n(-  +  —)  (A-27) 

where  Anr  describes  the  nonradiative  generation  of  carriers  and  rnr  is  the 
nonradiative  lifetime.  These  nonradiative  rates  are  often  much  faster  than 
the  radiative  terms,  and  dominate  the  thermalization  of  excited  carriers  in 
the  detector.  In  this  case,  illumination  by  the  signal  produces  a  smaller 
increase  in  carrier  concentration  ns  and  the  signal  to  noise  ratio  is  reduced. 

We  can  use  the  principle  of  detailed  balance  to  find  the  ratio  ns/nt  of 
the  signal  to  the  thermally-generated  carrier  concentration  including  nonra¬ 
diative  effects.  From  Equation  (A-27)  we  find: 

ns  =  Js@ar  ( Al  Nq  'jTtot  (A-28) 


where  rtot  is  the  total  lifetime  including  both  radiative  and  nonradiative 
processes: 

rtot  rT  Tnr 

In  the  present  context,  the  principle  of  detailed  balance  states  that  the  radia¬ 
tive  processes  and  the  nonradiative  processes  in  Equation  (A-27)  individually 
balance  in  thermal  equilibrium.  Thus  the  rates  of  radiative  generation  and 
recombination  are  unchanged  by  the  addition  of  the  nonradiative  process, 
and  Equation  (A-23)  for  the  thermal  carrier  concentration  remains  valid. 


Thus  for  the  general  case  including  both  radiative  and  nonradiative 
thermalization  of  carriers,  the  ratio  ns/nt  of  the  signal  to  the  thermally- 
generated  carrier  concentration  is: 

=  iJ* 

nt 


ns  f  J s\  f 

nt  \JtJ  \rTJ 


(A-30) 
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where  r0/rr  is  the  quantum  efficiency  for  the  radiative  recombination  of  ex¬ 
cited  carriers,  and  Jt  is  the  flux  of  blackbody  radiation  from  a  source  at  the 
same  temperature  as  the  detector.  This  expression  is  very  useful,  because  it 
allows  us  to  compare  different  approaches  to  detector  design.  Equation  (A- 
30)  tells  us  that  the  highest  signal-to-noise  ratio  at  a  given  temperature,  and 
thus  the  highest  operating  temperatures,  are  obtained  from  detectors  with 
the  highest  radiative  efficiency.  Designs  such  as  extrinsic  photoconductors 
and  quantum  well  detectors,  for  which  the  capture  of  excited  carriers  can 
proceed  via  phonon  emission,  often  require  lower  operating  temperatures  to 
compensate  for  greater  nonradiative  rates.  The  absolute  value  of  the  lifetime 
does  not  directly  influence  the  operating  temperature.  In  fact  shorter  radia¬ 
tive  lifetimes  are  desirable,  because  they  correspond  to  stronger  absorption 
and  permit  the  construction  of  smaller  pixels.  Kinch  and  Yariv  have  used  a 
similar  approach  to  study  the  effect  of  lifetime  on  operating  temperature  for 
quantum  well  detectors  as  discussed  below. 

Entropy  Ionization 


If  the  density  of  excited  states  in  an  infrared  photodetector  is  much 
larger  than  the  density  of  ground  states,  the  ground  states  thermally  depop¬ 
ulate  at  a  lower  temperature  than  one  might  expect  from  a  simple  Boltzmann 
factor  argument.  This  “entropy  ionization”  can  have  an  important  influence 
on  the  operating  temperature  of  quantum  detectors.  For  our  simple  two-level 
model,  the  volume  density  of  states  N2  for  the  upper  level  corresponds  to  the 
effective  density  of  states  in  the  conduction  band  N2  =  Nc: 

3/2 


Nc  =  2  n 


2i rm, 
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(A-31) 


where  is  the  number  of  conduction  band  minima,  and  me  is  the  electron 
density  of  states  effective  mass.  The  effective  density  of  states  can  be  quite 
large;  for  example,  with  nmin  =  l,me  =  0.1mo  and  T  =  100K,  we  find 
Nc  =  1.5  x  1017  cm-3.  For  extrinsic  photoconductors,  the  lower  level  in  our 
model  corresponds  to  donors  states  with  density  N  =  No,  and  the  traps 
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are  compensating  acceptors  with  density  No  =  N a-  Depending  on  the  type 
and  design  of  the  detector,  the  relative  density  of  levels  No,  Ni,  and  N2  vary 
by  orders  of  magnitude,  and  have  a  large  influence  on  the  required  operating 
temperature. 

At  low  temperatures  the  occupation  of  the  upper  level  is  small;  the 
thermal  carrier  density  is  nondegenerate  and  given  by  Maxwell-Boltzmann 
statistics: 

nt  =  N2e-{E2~^/kBT  (A-32) 

where  /j,  is  the  chemical  potential  (Fermi  level).  The  density  of  electrons  rq 
in  the  lower  level  is  degenerate  and  given  by  Fermi-Dirac  statistics: 

711  =  Nl  e{Ei-n)IkBT  +  i  (A-ZZ) 

The  total  carrier  concentration  is  fixed  and  equal  to  the  density  Ni  of  lower 
levels  minus  the  density  No  of  traps.  Setting  nt  +  n\  =  N\  —  N0  we  obtain  a 
quadratic  equation  for  the  chemical  potential,  and  thus  the  thermal  carrier 
concentration  nt.  At  temperatures  sufficiently  low  that  the  carrier  concen¬ 
tration  is  small  compared  with  the  trap  density  nt  <C  No  the  solution  of  this 
quadratic  equation  yields  approximately: 

nt  =  (Nt  -  No)  (j^j  e~AE'k*T  (A-34) 

the  same  as  Equation  (A-24)  above.  Thus  the  density  of  thermally-excited 
carriers  nt  is  equal  to  the  total  electron  density  (N\  —  No)  multiplied  by  a 
Boltzmann  factor  containing  the  energy  difference  A E  —  E2  —  E\  between 
the  upper  and  lower  levels,  weighted  by  the  relative  densities  N2jNo  of  the 
upper  level  and  traps. 

It  is  convenient  to  write  Equation  (A-34)  in  terms  of  the  change  in 
Helmholtz  free  energy  per  electron  A F  =  A E  —  TAS  between  the  upper 
and  lower  levels,  where  AS  is  the  difference  in  entropy  per  electron  between 
the  upper  and  lower  levels.  The  entropy  per  particle  is  approximately  the 
logarithm  of  the  number  of  available  states:  N2  for  the  upper  level  which  is 
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neaxly  empty,  and  7V0  for  the  lower  level  which  is  partially  full,  and  we  have 
approximately: 

AF  =  AE  -  kBTln  .  (A-35) 

Using  the  free  energy  difference,  the  thermally-excited  carrier  density  takes 
the  simple  form: 

nt  =  (M  -  N0)e~AF/kBT.  (A-36) 

That  is:  the  thermally-excited  carrier  concentration  is  just  the  product  of 
the  total  concentration  of  electrons  and  a  Boltzmann  factor  with  the  change 
in  free  energy  A F  taking  the  place  of  the  change  in  energy  AE.  This  use 
of  the  Helmholtz  free  energy  to  represent  the  internal  degrees  of  freedom  is 
standard  in  statistical  physics,  and  very  convenient  when  treating  the  relative 
occupation  of  levels  with  different  degeneracies. 


If  the  density  of  upper  states  is  larger  than  the  density  of  lower  states, 
carriers  gain  entropy  on  excitation,  and  the  number  of  excited  carriers  can  be 
much  larger  than  one  might  expect.  Furthermore,  the  lower  level  will  depop¬ 
ulate  at  a  lower  temperature  than  predicted  by  a  simple  Boltzmann  factor 
argument  using  the  energy  difference  between  levels  AE.  For  an  extrinsic 
photoconductor,  this  depopulation  corresponds  to  ionization  of  the  donor  lev¬ 
els.  We  can  estimate  the  ionization  temperature  7-  by  setting  A F/kBT  =  1 
in  Equation  (A-36): 


ksTi  = 


AE 

i+in(^y 


(A-37) 


Ordinarily  N2  is  a  function  of  temperature,  so  that  Equation  (A-37)  is  actu¬ 
ally  a  transcendental  equation  for  T{.  If  N2/N0  >  1,  as  is  often  the  case,  we 
find  that  the  lower  levels  ionize  at  temperatures  much  lower  than  those  corre¬ 
sponding  to  the  energy  difference  AE.  For  an  LWIR  detector  with  AE  =  100 
meV,  N2  =  1017  cm  3,  and  Nq  =  10l0  cm-3,  the  ionization  temperature  from 
Equation  (A-19)  is  T  ~  180K,  even  though  the  energy  gap  corresponds  to 
temperatures  T  ~  1000K.  Once  the  carriers  ionize,  the  free  carrier  concen¬ 
tration  becomes  temperature  independent,  and  the  device ^ is  insensitive  to 
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radiation.  In  order  to  operate  the  detector  successfully  it  must  be  cooled 
well  below  the  ionization  temperature. 

The  ionization  temperature  of  extrinsic  photoconductors  can  be  in¬ 
creased  by  increasing  the  trap  density  N0  to  a  significant  fraction  of  the 
density  of  lower  levels  iV2.  Increasing  N0  decreases  the  thermally-excited  car¬ 
rier  density  nth  (Equation  (A-34))  at  a  given  temperature  in  the  freeze-out 
regime  and  increases  the  ionization  temperature  2]  (Equation  (A-37)).  The 
largest  increase  in  ionization  temperature  is  obtained  when  Nq  approaches 
Ni,  and  a  significant  fraction  of  lower  states  are  empty  at  low  temperatures. 
In  general,  one  can  increase  the  ionization  temperature  by  making  Nq,  Ni, 
and  N2  all  the  same  order  of  magnitude,  thus  reducing  the  entropy  difference 
between  the  upper  and  lower  levels.  The  effect  of  trap  density  Nq  on  the  fig¬ 
ure  of  merit  ns/nt  for  noise  performance  (Equation  (A-30))  in  the  operating 
regime  below  the  ionization  temperature  is  determined  by  carrier  recombi¬ 
nation.  Increasing  Nq  proportionately  decreases  the  radiative  lifetime  rr  and 
can  also  change  the  nonradiative  recombination  rate.  Normally  a  reduction 
in  lifetime  is  undesirable,  because  the  change  in  carrier  density  ns  due  to 
the  signal  is  reduced.  However,  if  the  radiative  efficiency  is  high,  the  signal 
and  thermal  contributions  to  the  carrier  concentration  are  reduced  by  the 
same  factor,  and  the  ratio  ns/nt  remains  approximately  unchanged.  If  non¬ 
radiative  recombination  is  important,  then  the  dependence  of  the  radiative 
efficiency  on  trap  density  Nq  determines  noise  performance.  The  discussion 
above  assumes  that  the  detector  is  the  primary  source  of  noise  in  the  sys¬ 
tem.  If  the  detector  system  is  amplifier-noise  limited,  then  the  reduction  in 
photoconductive  gain  associated  with  a  shorter  carrier  lifetime  worsens  noise 
performance. 
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B  APPENDIX:  IMAGING  BOLOMETER 
ARRAYS 


Bolometers  operate  by  sensing  the  temperature  rise  associated  with  the 
absorption  of  radiation.  Bolometers  are  one  of  the  oldest  types  of  radiation 
detector,  and  they  have  many  advantages.  Bolometers  respond  to  absorbed 
energy,  and  can  be  made  sensitive  to  a  very  wide  range  of  wavelengths.  Be¬ 
cause  they  sense  heat  rather  than  photocarriers,  bolometers  are  insensitive 
to  the  photocarrier  population  and  dynamics.  As  a  result,  bolometers  can 
be  made  from  low-tech  materials,  and  they  can  be  operated  at  relatively 
high  temperatures.  These  properties  make  bolometers  well  suited  to  preci¬ 
sion  measurements  of  radiant  energy.  However,  conventional  bolometers  are 
typically  slow  and  insensitive  compared  with  quantum  detectors. 

Imaging  arrays  of  small  microbolometers  have  been  under  development 
since  the  1980s,  notably  through  the  HIDAD  (High  Density  Array  Devel¬ 
opment)  program  by  Texas  Instruments  and  by  Honeywell.  Using  mod¬ 
ern  semiconductor  processing  techniques,  large  staring  arrays  (238x340  pix¬ 
els)  of  small  (50x50  /j, m2)  bolometer  elements  were  fabricated  above,  and 
thermally  isolated  from,  a  Si  substrate  containing  readout  electronics.  The 
temperature-sensing  element  used  by  Texas  Instruments  is  a  ferroelectric  ma¬ 
terial,  while  that  used  in  the  Honeywell  approach  is  a  resistive  thermometer. 
In  the  early  1990s  the  program  was  declassified.  Since  then  development  of 
uncooled  microbolometer  arrays  has  been  actively  pursued  by  a  number  of 
companies  for  military  and  for  commercial  applications.  Santa  Barbara  Re¬ 
search  Center,  Lockheed  Martin,  Infrared  Solutions,  and  Raytheon-Amber 
are  using  variants  of  the  Honeywell  approach,  while  Texas  Instruments  has 
developed  a  second  generation  ferroelectric  approach.  These  bolometer  ar¬ 
rays  are  optimized  for  thermal  radiation  in  the  long-wave  infrared  (LWIR) 
atmospheric  transmission  band.  Good  noise  performance  at  room  temper- 
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ature  has  been  achieved  from  both  approaches.  Typical  values  of  the  noise 
effective  temperature  difference  in  1997  NETD  ~  40  mK  to  80  mK  for  f/1 
optics  at  video  frame  rates  (30  Hz  to  60  Hz).  The  ability  of  microbolometer 
arrays  to  operate  at  room  temperature  is  a  strong  advantage  over  quantum 
detectors,  which  must  be  cooled  to  ~  80  K.  However,  quantum  detector  ar¬ 
rays  made  from  HgCdTe  and  InSb  achieve  better  noise  performance  NETD 
~  10  mK  to  20  mK  in  the  mid-wave  (MWIR)  and  long-wave  infrared  atmo¬ 
spheric  transmission  bands. 

Microbolometer  arrays  will  likely  have  a  major  impact  on  infrared  imag¬ 
ing.  Their  sensitivity  and  image  quality  is  currently  adequate  for  many  appli¬ 
cations,  and  their  cost  could  fall  to  a  small  fraction  of  current  LWIR  imager 
prices  with  quantity  production.  Bolometers  are  inherently  multi-use  devices 
and  could  be  adapted  for  use  at  different  wavelengths  by  changing  the  ab¬ 
sorptive  coating.  The  performance  of  current  microbolometer  arrays  is  not 
near  their  fundamental  limits,  and  improvements  can  be  expected  with  con¬ 
tinued  development.  In  this  section  we  discuss  factors  which  currently  limit 
performance  and  make  suggestions  for  improvements. 

B.l  Bolometer  Noise 


A  bolometer  consists  of  a  temperature  sensing  element  with  heat  ca¬ 
pacity  C,  connected  to  a  heat  sink  through  a  thermal  conductance  G.  The 


device  senses  absorbed  radiation  via  the  increase  in  temperature  resulting 
from  the  deposited  energy.  The  temperature  rise  for  steady  illumination  is 
given  by: 


(B-l) 


where  P7  is  the  absorbed  power.  The  thermal  time  constant  of  the  bolometer 
is: 


r  = 


C 

G' 


(B-2) 
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The  bolometer  responds  to  changes  in  illumination  up  to  the  characteristic 
frequency  fc  =  l/2i rr,  above  fc  the  response  rolls  off. 

The  temperature  rise  AT  produced  by  illumination  is  typically  sensed 
via  a  corresponding  change  in  electrical  resistance  or  electric  polarization  of 
the  bolometer  element,  and  the  overall  response  is  proportional  to  the  tem¬ 
perature  coefficient  of  resistance  dR/dT  or  polarization.  High  temperature 
coefficients  are  necessary  to  overcome  electronic  noise  in  the  readout  electron¬ 
ics  and  the  bolometer  itself.  The  responsitivity  of  bolometers  is  discussed  in 
detail  below.  • 


The  fundamental  noise  limit  for  bolometers  is  set  by  temperature  fluc¬ 
tuations  in  the  sensitive  element,  given  by: 

(ST2)  =  (B-3) 

where  T  is  the  temperature  of  the  bolometer.  For  a  given  deposited  energy 
A E,  the  ratio  of  the  rms  thermal  fluctuation  to  the  temperature  rise  is: 

Thus  the  best  energy  resolution  is  obtained  with  the  smallest  heat  capacity. 
Small  heat  capacity  also  reduces  the  bolometer  time  constant  and  improves 
the  response  speed. 


It  is  interesting  to  relate  thermal  fluctuation  noise  to  the  number  of 
atoms  Na  in  the  sensitive  element.  Above  the  Debye  temperature  Tp  of  the 
bolometer  material,  all  of  the  thermal  degrees  of  freedom  are  excited,  and  the 
heat  capacity  is  C  =  3 Naks  from  the  principle  of  equipartition.  For  this  case 
bolometer  noise  is  proportional  to  the  square  root  of  the  number  of  degrees 
of  freedom: 


(6T2)1/2 

AT 


kBT 
A  E 


(3  Na)1'2. 


(B-5) 


To  reduce  temperature  fluctuations  above  the  Debye  temperature  one  must 
reduce  the  number  of  atoms  in  the  bolometer  element  or  reduce  the  tem¬ 
perature.  Below  the  Debye  temperature  the  heat  capacity  decreases  rapidly 
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as  vibrational  motion  is  frozen  out,  and  noise  performance  improves  corre¬ 
spondingly.  Thus  it  is  desirable  to  choose  bolometer  materials  with  Debye 
temperatures  larger  than  the  temperature  of  operation  Tp  >  T .  if  possible. 
For  Si  the  Debye  temperature  is  7p  =  645  K. 

The  power  spectrum  of  thermal  fluctuations  in  the  bolometer  element 
can  be  found  from  the  autocorrelation  function,  which  decays  exponentially 
in  time  with  the  bolometer  time  constant  r: 


(6T(t0)ST(t0  + 1))  =  (6T2)exp(—\t\/r).  (B-6) 


The  power  spectrum  for  temperature  fluctuations  is  proportional  to  the 
Fourier  transform  of  the  correlation  function: 


/OO 

(5T(t0)5T(t0  +  t))e~^tdt. 

-oo 


Performing  the  integral  we  have: 


1 

1  +  (27 r/r)2 ' 


(B-7) 


(B-8) 


Below  the  characteristic  frequency  fc  =  1/27 rr,  the  power  spectrum  is 
frequency  independent,  and  the  noise  equivalent  illumination  power  due  to 
thermal  fluctuations  is  given  by: 


<«a> - 


(B-9) 


where  rj  is  the  quantum  efficiency  for  absorption  of  photons.  From  Equation 
(B-9)  we  find  that  the  noise  power  can  be  reduced  by  reducing  the  ther¬ 
mal  conductance  G  and  by  reducing  the  temperature.  The  heat  capacity 
does  not  appear  in  this  expression,  because  the  temperature  fluctuations  are 
spread  over  a  bandwidth  proportional  to  1/C.  Combining  these  results,  we 
find  that  bolometer  noise  performance  improves  as  we  reduce  the  heat  ca¬ 
pacity  and  thermal  conductance.  Thus  bolometers  are  natural  candidates  for 
microfabrication. 
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B.2  Comparison  of  Bolometers  with  Quantum  Detec¬ 
tors 


Photon  shot  noise  determines  the  fundamental  background-limited  in¬ 
frared  performance  (BLIP)  for  any  type  of  infrared  detector.  The  average 
fluctuation  8P~f  in  illumination  power  incident  on  the  detector  due  to  photon 
shot  noise  is: 

(6P72)  =  2PyhuAf  (B-10) 


where  P7  is  the  illumination  power  and  hu  is  the  photon  energy.  The  ratio 
of  shot  noise  to  thermal  fluctuation  noise  in  a  bolometer  is: 


(J>Pt)  =  i  [A r7\  (  hv  \ 
(5P?h)  2  V  T  )  \kBT ) 


(B-ll) 


where  A T7  =  'qP^/G  is  the  temperature  increase  produced  by  illumination. 
This  relation  shows  that  it  is  difficult  to  make  infrared  bolometers  which  ap¬ 
proach  background  limited  performance  for  low  illumination  levels  for  which 
ATy  is  small.  Bolometers  are  well  suited  to  precise  measurements  of  large 
signals,  such  as  the  thermal  emission  of  300  K  blackbodies,  and  can  ap¬ 
proach  the  BLIP  limit  in  this  case.  For  uncooled  microbolometers  looking 
at  300  K  sources  in  the  LWIR  band,  hv/kBT  4,  A  r7  1  K  (see  below), 
and  (SP])/(6PS,)  ~  0.01,  so  that  the  theoretically  detectable  temperature 
difference  NETD  ~  10  mK  is  ~  10  times  larger  than  the  BLIP  limit.  Cool¬ 
ing  and  operation  in  shorter  wavelength  bands  improves  performance  pro¬ 
portionately.  For  visible  light  hu/kBT  ~  80  (T  =  300  K),  and  uncooled 
bolometer  arrays  are  theoretically  capable  of  approaching  background  lim¬ 
ited  performance  (JP7) /  (5P?h)  ~  1  for  high  levels  of  illumination  such  as 
daytime  images  of  the  earth’s  surface. 


The  sensitivity  and  noise  performance  of  bolometers  improves  as  the 
thermal  conductance  G  decreases.  A  lower  limit  to  G  is  determined  by  ther¬ 
mal  coupling  to  the  detector  environment  via  blackbody  radiation.  From  the 
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Stefan-Boltzmann  law  we  find: 


Gbb  =  4£(7sbT3  A 


(B-12) 


where  s  and  A  are  the  emissivity  and  area  of  the  bolometer  element,  and 
&sb  —  n2kB/6Qh3c2  =  5.67  x  10-8  W/m2K4  is  the  Stefan-Boltzmann  con¬ 
stant.  For  HIDAD  specifications  we  find  Gbb  —  1.5  x  10-8  W/K.  The  tem¬ 
perature  rise  AT7  produced  by  illumination  in  the  radiative  coupling  limit 
is: 


A  T7 
A  Ts 


~  I  (-)  (±)  (h\ 

4  v 7r  y  Us/  V  t  ) 


(B-13) 


where  Q,  is  the  solid  angle  through  which  the  detector  is  illuminated  by  a 
source  of  emissivity  es  at  temperature  Ts .  For  uncooled  microbolometer  ar¬ 
rays  the  source  and  detector  temperatures  are  the  same  and  ATy/ATs  = 
D/ 4tt,  giving  A T7  ~  20  K  for  f/1  optics,  which  corresponds  to  a  noise  ef¬ 
fective  temperature  difference  NETD  ~  4  mK.  Cooling  reduces  the  thermal 
conductance  due  to  radiative  coupling  dramatically.  For  example,  cooling 
from  room  temperature  to  80  K  reduces  Gbb  by  a  factor  ~  50. 


We  can  compare  the  performance  of  radiatively  coupled  bolometers  to 
quantum  detectors  for  wide  band  blackbody  sources  by  inserting  the  temper¬ 
ature  rise  from  the  previous  expression  into  the  ratio  of  shot  to  temperature 
fluctuation  noise  powers: 


{^)  _  1  (S}\  fes\  ( Ts\ 4  (  hu  \ 

(5  Pth)  8WUA  t)  \kBT) 


(B-14) 


where  hv  is  the  average  photon  energy  for  the  blackbody  source.  Doing  the 
average  for  shot  noise  distributed  over  the  blackbody  spectrum,  and  assuming 
nondegenerate  statistics,  one  finds: 


hu  = 


360C(5) 


kBT  ~  3.83 kBT 


(B-15) 


where  ^(5)  is  the  Riemann  zeta  function  for  argument  5. 


Radiative  coupling  contributes  an  additional  source  of  noise,  which  is 
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simply  the  shot  noise  associated  with  the  absorption  and  emission  of  black- 
body  photons  at  the  detector  temperature  T: 

($P‘bb)  =  4Pj5ghi/A/  (B-16) 


where  Pbb  =  AcosbT 4  is  the  total  power  radiated  by  the  bolometer  element, 
and  hv  is  given  above.  The  numerical  factor  is  4  instead  of  2,  because  both 
absorption  and  emission  of  photons  contribute  to  the  noise.  We  can  assess 
the  importance  of  this  noise  source  for  detectors  in  which  blackbody  radia¬ 
tion  determines  the  thermal  conductance  by  taking  the  ratio  of  shot  noise 
associated  with  the  signal  to  temperature  fluctuation  noise  in  the  radiation 
coupled  limit: 


(«?)  _  1  (Sl\  /£S\  (Ts\* 

(SPbb)  2  \7T  /  \  £  )  V  T  ) 


(B-17) 


For  equal  emissivities  and  temperatures  in  uncooled  detectors  this  ratio  is 
(SPy) / (SPqB)  =  Q/27T,  always  less  than  unity  (~1/10  for  f/1  optics),  be¬ 
cause  the  solid  angle  associated  with  the  source  is  smaller  than  27t.  Thus  for 
radiatively  coupled  uncooled  bolometers  the  noise  performance  cannot  reach 
the  BLIP  limit.  If  the  detector  is  cooled,  the  shot  noise  associated  with 
radiation  coupling  falls  rapidly,  and  the  BLIP  limit  is  reached  at  moderate 
temperatures  (T  ~  200  K  for  f/1  optics).  Whenever  the  thermal  conductance 
G  exceeds  the  radiation-coupled  limit,  conventional  thermal  fluctuation  noise 
dominates. 


B.3  Responsitivity 


The  noise  equivalent  power  5P1  of  an  ideal  bolometer  is  given  by: 

(SPy)  =  4kBT2GAf  (B-18) 

where  T  is  the  temperature,  G  is  the  thermal  conductance  to  the  substrate, 
and  A /  is  the  bandwidth  in  Hz.  The  response  A R  of  a  resistive  bolometer  is 
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proportional  to  the  temperature  coefficient  dR/dT,  which  can  be  expressed 
in  dimensionless  form: 


dlogR  _  /T\  dR 
dlogT  ~  vpj  dT' 


(B-19) 


the  temperature  coefficient  a  is  the  slope  of  a  log-log  plot  of  R  vs.  T.  Then 
the  resistance  change  corresponding  to  the  absorbed  power  P7  is: 


A  R  __  a  P7 
~R  =  ~GT' 


(B-20) 


The  resistance  fluctuations  SR  due  to  temperature  fluctuations  in  the  bolome¬ 
ter  are  given  by: 

(fff) <*»> 

For  temperature  fluctuation  noise  to  dominate  in  a  resistive  bolometer,  it 
must  be  larger  than  Johnson  noise.  Assuming  that  a  bias  current  IB  is  used 
to  measure  the  resistance,  the  ratio  of  these  two  noise  sources  is  given  by: 


4tfI2BR2kBAf/G  o?I%R  2/A  Tb\ 

(SVJ)  iRkBTAf  GT  a  {  T  ) 


where  A TB  is  the  heating  associated  with  the  bias  current  IB.  For  this  ratio 
to  be  greater  than  one,  the  dimensionless  temperature  coefficient  should  be 
large  a  »  1. 


Using  this  notation,  we  can  compare  several  candidate  materials  for  re¬ 
sistance  bolometers.  Simple  metals  at  high  temperatures  show  approximately 
power  law  temperature  dependence: 


R  oc  Ta.  (B-23) 

For  this  case  the  temperature  coefficient  is  simply: 

a  =  a.  (B-24) 

It  is  difficult  to  make  sensitive  bolometers  from  simple  metals,  because  the 
exponent  a  is  not  large:  typically  a  ~  1. 
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The  resistance  of  a  common  class  of  bolometer  materials  including  car¬ 
bon  and  vanadium  oxide  is  characterized  by  variable  range  hopping  for  which: 


R  (x  exp  [C T0/T )1/4] .  (B-25) 


For  this  case  the  dimensionless  temperature  coefficient  is: 

1/4 


(B-26) 


For  variable  range  hopping  it  is  difficult  to  achieve  large  temperature  coeffi¬ 
cients  a  even  if  T0/T  is  large,  due  to  the  exponent  1/4. 


Semiconductors  are  another  candidate  material.  For  this  case  the  resis¬ 
tance  decreases  with  temperature  as: 


R  oc  exp {Eact/kBT) 


(B-27) 


where  Eact  is  the  activation  energy  for  free  carriers.  The  dimensionless  tem¬ 
perature  coefficient  is: 


a  = 


Jact 


kBT 


(B-28) 


For  extrinsic  semiconductors  at  low  temperatures  in  the  freezeout  regime  the 
activation  energy  is  the  donor  or  acceptor  binding  energy  Eact  =  EB-  For  this 
case  it  is  possible  to  achieve  very  large  dimensionless  temperature  coefficients 
at  very  low  temperatures,  for  example  a  ~  100  for  acceptors  in  Ge  at  T  = 
0.3  K,  and  very  sensitive  bolometers  have  been  demonstrated. 


For  intrinsic  semiconductors  operating  at  higher  temperatures,  the  ac¬ 
tivation  energy  is  half  the  energy  gap,  Eact  =  Egap/2.  A  lower  limit  to  the 
operating  temperature  and  an  upper  limit  to  a  in  this  case  is  given  by  the 
transition  to  extrinsic  behavior  when  the  intrinsic  carrier  concentration  be¬ 
comes  comparable  to  the  impurity  concentration.  For  intrinsic  Ge  at  T  = 
300  K,  the  dimensionless  temperature  coefficient  a  ~  10. 


A  large  temperature  coefficient  can  be  obtained  at  a  particular  tem¬ 
perature  by  selecting  a  bolometer  material  which  undergoes  an  electronic 
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phase  transition.  High  temperature  superconductors  are  being  developed 
as  bolometers,  because  the  resistance  rapidly  drops  near  the  superconduct¬ 
ing  transition  temperature  Tc.  For  YBCO,  the  transition  temperatures  Tc 
can  be  made  comparable  to  standard  cooler  temperatures  T  80  K,  and 
temperature  coefficients  higher  than  conventional  bolometer  materials  are 
possible.  Honeywell  has  reported  high  Tc  bolometers  with  temperature  co¬ 
efficients  (1  /R)dR/dT  =  0.2  K"1  near  T  =  70  K,  for  which  ct  ~  14.  The 
resistive  transition  in  this  material  is  broad  in  order  to  reduce  the  require¬ 
ment  for  temperature  regulation.  Much  sharper  transitions  are  obtained  in 
high  quality  samples  of  high  Tc  materials.  It  is  not  necessary  to  use  super¬ 
conductors  for  bolometers,  because  only  the  resistance  change  is  important. 
Other  resistive  transitions  can  be  exploited.  For  example,  the  electrons  in 
NbSe3  undergo  a  phase  transition  to  a  charge  density  wave  state  at  moderate 
temperatures,  accompanied  by  a  sharp  increase  in  resistance. 
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C  APPENDIX:  THERMOELECTRIC 
MATERIALS 


Thermoelectric  materials  are  those  in  which  the  flow  of  heat  is  coupled 
to  the  electric  current.  Because  electrons  carry  entropy  as  well  as  electrical 
charge,  all  normal  conductors  are  thermoelectric  to  some  extent.  Good  ther¬ 
moelectric  materials  with  relatively  strong  coupling  between  the  thermal  and 
electric  currents  are  useful  for  applications  ranging  from  solid  state  coolers, 
to  temperature  sensors,  to  power  generation. 

For  this  report,  the  most  relevant  applications  of  thermoelectric  ma¬ 
terials  are  in  solid  state  coolers  and  bolometric  infrared  detectors.  Solid 
state  coolers  using  thermoelectric  elements  are  compact  and  have  no  moving 
parts.  However,  their  efficiency  is  relatively  low,  and  present  materials  reach 
temperatures  only  ~  200  K.  Thermoelectric  materials  with  better  low  tem¬ 
perature  performance  are  needed  to  push  the  operating  temperature  lower. 
Thermoelectric  materials  can  also  be  used  for  bolometric  infrared  sensors. 
The  best  thermoelectric  materials  perform  better  than  simple  thermocouples 
for  this  application,  and  uncooled  microbolometer  arrays  using  this  approach 
are  under  development  at  the  Jet  Propulsion  Laboratory.  High  temperature 
thermoelectrics  operating  at  temperatures  ~  800  K  can  be  used  to  generate 
electric  power. 

New  low  temperature  thermoelectric  materials  are  needed  to  make  cool¬ 
ers  capable  of  reaching  temperatures  below  200  K.  Because  the  characteristics 
of  thermoelectric  materials  have  been  thoroughly  explored  since  the  1950s, 
this  task  is  difficult.  However  the  development  of  new  classes  of  man-made 
materials  -  heterostructures  and  nanostructured  materials  -  creates  hope 
that  better  low  temperature  thermoelectrics  may  be  discovered  or  designed. 
In  order  to  provide  a  guide  to  desirable  physical  properties  for  new  thermo¬ 
electric  materials,  we  review  the  origins  of  thermoelectricity  below. 
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To  understand  the  physical  mechanism  responsible  for  thermoelectricity 
and  to  compare  different  materials,  it  is  helpful  to  develop  a  simple  model. 
Practical  thermoelectrics  are  typically  semimetals  or  heavily  doped  semi¬ 
conductors.  We  consider  a  material  with  one  type  of  carrier  (electrons  or 
holes)  with  charge  q,  effective  mass  m,  and  cubic  crystal  symmetry  so  that 
the  electrical  and  thermal  conductivity  are  isotropic.  Detailed  treatments  of 
thermoelectricity  including  specific  material  properties  are  given  in  reviews 
(Mahan  1997,  CRC  Handbook  of  Thermoelectrics  1995)  and  textbooks  (see 
Ziman  1972). 


The  thermoelectric  effect  occurs  in  all  normal  conductors  because  elec¬ 
trons  (or  holes)  carry  both  electrical  charge  and  heat.  The  entropy  s  per 
charge  carrier  is  quite  generally  given  by  the  expression: 


where  (e)  is  the  energy  per  particle,  p  is  the  chemical  potential,  and  T 
is  the  absolute  temperature.  Equation  (C-l)  applies  equally  well  to  low 
density  (non-degenerate)  carrier  gases  in  semiconductors  and  high  density 
(degenerate)  gases  in  metals.  This  expression  arises  in  the  derivation  of  the 
Boltzmann  factor  and  can  be  found  from  the  Helmholtz  free  energy  : 


s  = 


1_ 

T 


(g)  -  At 
T 


(C-2) 


Here  U  and  Ss  are  the  total  thermal  energy  and  entropy  of  a  gas  of  N  particles 
in  volume  V. 


Because  thermoelectricity  is  a  transport  phenomenon  we  consider  the 
nonequilibrium  distribution  f(p)  of  an  ideal  gas  of  particles  under  the  ac¬ 
tion  of  an  applied  electric  field  E,  a  gradient  in  chemical  potential  V/r,  and 
a  temperature  gradient  VT.  In  the  relaxation  time  approximation  of  the 
Boltzmann  transport  equation  (see  Ziman  1972): 
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Here  e  and  v  are  the  particle  energy  and  velocity  for  momentum  p.  The 
three  terms  inside  the  square  brackets  are  the  three  forces  which  produce 
particle  transport:  the  electric  force  qE,  the  diffusive  force  due  to  density 
gradients  —Vp,  and  a  thermal  force  produced  by  temperature  gradients 
— [(e  —  p)/T)VT.  It  is  the  thermal  force  which  is  responsible  for  thermo¬ 
electricity. 


Equation  (C-3)  assumes  that  the  nonequilibrium  distribution  f(p)  of 
particles  in  momentum  space  remains  near  the  thermal  distribution  f°(p) 
under  the  action  of  the  applied  forces,  so  that  the  change  can  be  expressed 
as  f—f°  =  (—df°/de)Ae.  The  energy  Ae  imparted  to  carriers  of  momentum 
p  by  the  force  is  the  product  of  the  force,  the  carrier  velocity,  and  the  mo¬ 
mentum  relaxation  time  rp  as  in  Equation  (C-3).  We  have  also  assumed  that 
the  collisions  which  relax  the  distribution  are  elastic,  or  energy  conserving. 
An  example  of  an  elastic  process  is  impurity  scattering;  examples  of  inelas¬ 
tic  processes  are  electron-electron  and  electron-phonon  scattering.  Equation 
(C-3)  can  be  used  as  the  starting  point  for  treating  carrier  transport  in  both 
metals  and  semiconductors.  However,  good  metals  generally  make  poor  ther¬ 
moelectric  materials,  so  we  assume  that  the  particle  gas  is  non-degenerate, 
i.e.  that  is  a  low  density  gas  obeying  Maxwell  Boltzmann  statistics: 

/"  =  «pexp  (^)  with  ?£  =  -jL;  f  (CM) 

where  Np  is  the  density  of  states  in  p-space. 


One  can  understand  the  origin  of  the  thermal  force  by  noting  that  par¬ 
ticles  move  more  rapidly  in  hotter  regions  of  the  sample,  and  thus  tend  to 
escape  a  given  region  faster  than  they  are  replaced  by  cold  particles,  pro¬ 
ducing  a  particle  current  flowing  from  hot  to  cold.  The  coefficient  which 
determines  the  strength  of  thermal  force  is  simply  the  entropy  per  particle 
s.  For  a  non-degenerate  gas: 


s  = 


£  —  p 

T  ‘ 


(C-5) 
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It  is  helpful  to  note  that  s  is  quite  generally  related  to  the  chemical  potential 
via  the  thermodynamic  identity: 

-(»)„-(*).  ~(*L-  m 

Here  Ss  is  the  total  entropy  for  N  particles  in  volume  V,  and  F  is  the 
Helmholtz  free  energy. 


The  electric  current  density  J  produced  by  the  sum  of  these  forces  is 
determined  by  the  integral  over  momentum  space: 

2  q 
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One  can  gain  insight  into  the  physical  origins  of  thermoelectricity  without 
actually  carrying  out  the  integration.  For  isotropic  materials  the  total  electric 
current  has  the  following  form: 


J  =  a{E-  V//)  -  aSVT  (C-8) 

where  cr  is  the  electrical  conductivity,  and  S  is  the  Seebeck  coefficient,  a  ma¬ 
terial  parameter  which  characterizes  the  thermoelectric  current.  The  electric 
and  diffusive  forces  are  grouped  together,  because  it  is  generally  their  sum 
which  is  measured  in  experiments  (consider  for  example  the  ‘voltage’  across 
a  pn  junction  diode). 


The  electrical  conductivity  is  given  by 

2q2 


kBT 


fd3Pf° 


rpvv  •  e 


(C-9) 


where  e  is  a  unit  vector  in  the  direction  of  the  electric  field.  This  expression 
can  treat  complex  situations  with  momentum  dependent  relaxation  times  Tp 
and  densities  of  states  Np.  One  generally  side  steps  these  complications  by 
defining  a  single  effective  scattering  time  r  and  an  effective  mass  m,  which 
are  related  to  microscopic  processes  via  the  integral  above.  We  then  get  the 
well  known  expression: 

nq2T 

o’  = -  (C-10) 
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where  n  is  the  carrier  concentration,  and  we  have  used  {v£)  =  ksT/m  for  the 
average  component  of  velocity  along  the  electric  field. 

In  the  same  spirit,  we  can  obtain  an  expression  for  the  Seebeck  coefficient 
S.  The  integral  for  oS  is  identical  to  that  for  the  electrical  conductivity  a , 
except  for  the  additional  factor  (s  —  n)/T.  Assuming  that  the  integrand 
is  well-behaved,  we  simply  move  the  physically  relevant  quantities  outside, 
replacing  the  energy  by  the  effective  average  value  e,  and  find: 

HWK  <«■> 

As  shown,  the  Seebeck  coefficient  is  equal  to  the  entropy  per  free  carrier  $ 
divided  by  the  charge  q;  S  can  have  either  sign,  and  it  is  positive  for  positively 
charged  carriers.  This  is  a  simple  and  important  result  which  will  allow  us  to 
assess  the  relative  merits  of  different  thermoelectric  materials.  However,  one 
should  keep  in  mind  that  useful  thermoelectric  materials  may  rely  on  more 
complex  microscopic  processes  which  violate  simple  rules. 

A  simple  physical  consequence  of  the  thermal  force  can  be  found  by 
evaluating  Equation  (C-8)  for  the  situation  in  which  a  temperature  gradient 
is  set  by  a  flow  of  heat  through  the  material,  but  no  electric  current  can  flow: 

E  -  if/i  =  SVT.  (C-12) 

q 

We  see  that  thermal  gradients  in  the  absence  of  an  electric  current  produce 
electric  fields  (more  precisely  the  electrochemical  field  E  —  (l/q)V//);  the 
constant  of  proportionality  is  the  Seebeck  coefficient  S.  Hot  carriers  flow 
away  from  higher  temperature  regions  of  the  sample  more  rapidly  than  they 
are  replaced  by  colder  carriers.  This  flow  is  offset  by  the  electrochemical  field 
coming  from  the  build  up  of  charge  ( E )  and  carriers  (— V/u)  in  colder  parts 
of  the  sample.  If  the  carriers  cannot  flow  out  of  the  sample,  the  total  current 
density  must  cancel  in  equilibrium  leading  to  Equation  (C-12).  For  metals 
the  resulting  electric  field  E  is  typically  much  greater  than  the  diffusive 
current  and  one  can  neglect  the  diffusive  force  —  V/i  in  Equation  (C-12);  for 
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semiconductors  inclusion  of  the  diffusive  force  may  be  necessary  (e.g.  for  the 
thermoelectric  effect  in  inhomogeneous  devices  like  pn  junctions). 

Just  as  a  temperature  gradient  produces  an  electric  current,  an  electric 
field  produces  a  thermal  current.  For  linear  response  the  thermal  current 
density  Jq  can  be  written  in  the  form: 

JQ  =  aST  (e  -  iv/ij  -  KqEVT  (C-13) 

where  KQE  is  the  thermal  conductivity  for  E  =  0.  Note  that  KqE  is  not  the 
usual  thermal  conductivity  Kq  that  one  measures  in  the  absence  of  electric 
currents;  we  return  to  this  point  later.  The  coefficient  of  the  electrochemical 
field  contains  precisely  the  same  factor  aS  that  appears  in  Equation  (C-8) 
for  the  electrical  current  density.  This  is  no  accident,  but  an  exact  symmetry 
which  follows  from  time-reversal  symmetry  in  the  underlying  equations  of 
motion,  an  example  of  an  Onsager  relation  (see  Landau  and  Lifshitz  1969). 

One  can  understand  the  physical  origin  of  heat  transport  by  carriers  by 
noting  that  the  heat  carried  per  particle  in  a  flow  of  charge  carriers  is: 

A  Q  =  Ts  =  qST  (C-14) 

where  s  is  the  entropy  per  carrier,  and  we  have  used  Equation  (C-ll)  for  the 
Seebeck  coefficient  S.  Thus  the  heat  current  produced  by  an  electric  current 
with  VT  =  0  is  given  by: 

JQ  =  STJ  =  a  ST  {e  -  ^  (C-15) 

in  agreement  with  Equation  (C-13). 

Equations  (C-8)  and  (C-13)  together  are  a  starting  point  for  the  theory 
of  thermoelectric  devices  (see  Mahan  1997).  For  applications,  the  natural 
variables  are  J  and  VT,  rather  than  E  —  (l/^Vyu  and  VT  .  Eliminating 
E  -  (1  /q)Vfj,  from  Equation  (C-13)  using  Equation  (C-8)  we  find: 

Jq  =  STJ  -  KqVT  (C-16) 
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where  Kq  is  the  thermal  conductivity  for  zero  current  density  J  =  0: 

Kq  =  Kqe  1  -  (^)  T  .  (C-17) 

This  is  the  usual  thermal  conductivity  that  one  determines  in  purely  thermal 
measurements.  The  thermal  conductivity  Kq  is  smaller  than  Kqe,  because 
the  heat  carried  by  the  particle  current  is  subtracted. 

An  important  figure  of  merit  for  thermoelectric  materials  is  the  quantity 
Z  which  has  dimensions  of  inverse  temperature: 

<c-18> 

The  efficiency  of  thermoelectric  coolers  is  determined  by  the  dimensionless 
product  ZT  which  must  be  comparable  or  greater  than  one  for  good  per¬ 
formance  (see  Mahan  1997).  In  the  limit  ZT  »  1  the  cooling  efficiency 
approaches  the  ideal  Carnot  efficiency.  The  best  current  thermoelectric  ma¬ 
terials  have  ZT 


We  are  now  in  a  position  to  evaluate  the  promise  of  thermoelectric 
materials  based  on  their  physical  characteristics.  A  good  starting  point  is 
the  Seebeck  coefficient  for  nondegenerate  carriers: 


kBf  3 fJ_\ 

q  \  2  kBTJ 


(C-19) 


where  we  have  used  e  =  ( 3/2)kBT ,  and  the  constant  kB/q  —  86/uV/K.  The 


chemical  potential  is: 


/i  =  kBT\n 
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\  2ttH2  ) 


(C-20) 


Here  uq  is  the  quantum  concentration  (effective  density  of  states),  with  g 
the  product  of  the  spin  and  valley  degeneracies.  For  example  g  =  12  for 
n-type  Si,  which  has  spin  1/2  and  six  electron  valleys.  By  definition,  n  <  uq 
for  nondegenerate  carriers,  and  g,  <  0.  The  chemical  potential  for  semicon¬ 
ductors  is  determined  by  the  dopant  concentrations  as  well  as  the  quantum 
concentration  (see  Seeger  1988,  Sze  1981). 
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It  is  instructive  to  consider  the  temperature  dependence  of  the  Seebeck 
coefficient  through  the  chemical  potential  for  an  n-type  semiconductor  with 
donor  concentration  Njq.  At  high  temperatures  the  material  is  intrinsic  with 
equal  electron  and  hole  concentrations,  and  the  chemical  potential  lies  near 
the  middle  of  the  energy  gap  eg.  Intrinsic  semiconductors  are  not  typically 
good  thermoelectrics,  because  the  Seebeck  coefficients  of  electrons  and  holes 
subtract  and  tend  to  cancel.  As  the  temperature  falls  to  values  ksT  <  eg/2 
the  chemical  potential  rises  toward  the  conduction  band,  and  holes  become 
the  minority  carrier.  In  this  extrinsic  regime,  most  useful  for  thermoelectrics, 
the  electron  concentration  is  constant  and  equal  to  the  donor  concentration. 
At  sufficiently  low  temperatures  the  chemical  potential  becomes  pinned  at 
the  donor  binding  energy  —B,  and  the  carrier  concentration  freezes  out  ex¬ 
ponentially.  From  this  development  we  see  that  the  Seebeck  coefficient  for 
nondegenerately  doped  semiconductors  increases  with  temperature  until  the 
transition  to  intrinsic  behavior  is  approached.  Mahan  (1997)  gives  an  em¬ 
pirical  rule  that  the  best  operating  temperature  for  semiconductor  thermo¬ 
electrics  is  ksT  ~  £g/ 10  in  agreement  with  these  arguments. 


Good  thermoelectric  materials  must  have  high  electrical  conductivity  a 
as  well  as  a  large  Seebeck  coefficient  5,  so  that  the  figure  of  merit  ZT  can  be 
as  large  as  possible.  High  electrical  conductivity  increases  the  cooling  power 
by  increasing  the  flow  of  particles  (and  entropy)  for  a  given  applied  voltage. 
The  product  appearing  in  the  numerator  of  Z  is  called  the  power  factor: 


(C-21) 


Assuming  t  is  density  independent,  the  power  factor  has  a  maximum  at 
In (n/nQ)  =  -1/2  for  which  the  carrier  density  n  is  comparable  to  the  quan¬ 
tum  concentration  uq.  The  condition  n  ~  uq  also  marks  the  boundary  above 
which  the  gas  of  carriers  becomes  degenerate  and  obeys  Fermi-Dirac  statis¬ 
tics.  Materials  having  this  property  are  heavily  doped  semiconductors  and 
semimetals.  The  magnitude  of  the  Seebeck  coefficient  at  maximum  power 
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factor  S  =  2ks/q  =  170yuV/K  is  typical  of  good  thermoelectric  materials 
(Mahan  1997). 

The  influence  of  the  effective  mass  m  on  the  performance  of  thermo¬ 
electric  materials  can  be  examined  by  evaluating  the  power  factor  at  its 
maximum  at  n  =  n^e-1/2: 

(crS2) max  =  ra1/2T3/2.  (C-22) 

As  shown,  the  maximum  power  factor  increases  with  the  effective  mass  m 
and  temperature  T,  because  the  quantum  concentration  is  larger  and  greater 
carrier  concentrations  are  possible.  A  heavy  effective  mass,  or  high  density  of 
states  at  the  Fermi  surface,  is  generally  desirable  for  thermoelectric  materials. 
Unfortunately,  small  bandgap  semiconductor  materials  which  could  be  useful 
thermoelectrics  at  room  temperature  and  below  tend  to  have  light  effective 
masses. 

Thermal  conductivity  is  the  final  major  factor  influencing  the  perfor¬ 
mance  of  thermoelectric  materials.  Low  thermal  conductivity  improves  the 
efficiency  of  thermoelectric  coolers  by  reducing  unwanted  reverse  heat  flow. 
The  total  thermal  conductivity  Kq  is  the  sum  of  contributions  from  charge 
carriers  and  from  phonons.  Because  phonons  do  not  contribute  to  the  ther¬ 
moelectric  effect  or  the  cooling  power,  the  best  performance  would  be  ob¬ 
tained  if  it  were  possible  to  eliminate  heat  conduction  by  phonons,  and  free 
carriers  were  responsible  for  the  transport  of  both  charge  and  heat.  For  good 
metals  with  n  ~  1023  cm-3  this  is  the  case,  but  for  heavily  doped  semiconduc¬ 
tors  and  semimetals  with  much  lower  carrier  concentrations  n  ~  1019  cm-3, 
phonons  typically  dominate  thermal  conduction.  Nonetheless  it  is  interest¬ 
ing  to  examine  the  effects  of  heat  conduction  by  free  carriers.  The  thermal 
conductivity  of  a  nondegenerate  carrier  gas  is  (see  Kittel  1971): 

Kq  =  \cv\  =  ~k%T  (C-23) 

where  C  is  the  heat  capacity  per  unit  volume.  Because  carriers  transport 
both  charge  and  heat,  the  ratio  of  thermal  to  electrical  conductivity  is  a 
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constant: 


2 


Kq  _  (3/2)fc|T  _  3  (kB\ 

o  —  _  2  J  T'  <C-24> 

This  relation  is  known  as  the  Wiedemann-Franz  law.  [For  good  metals 
it  takes  a  slightly  different  form  due  to  the  difference  in  carrier  statistics: 


Kq/u  =  (TT2/3)(kB/q)2T.]  Inserting  Equation  (C-24)  into  Equation  (C-18) 
we  get  an  upper  limit  to  the  figure  of  merit: 


(£)’■  (-) 
Using  typical  values  for  the  entropy  per  particle  s/kB  ~  1  in  this  expression, 
one  finds  that  values  ZT  >  1  are  difficult  to  obtain,  in  agreement  with 
practical  experience. 


Given  that  phonons  dominate  the  thermal  conductivity  Kq  of  most 
thermoelectric  materials,  one  would  like  to  reduce  Kq  by  decreasing  the 
population  of  phonons  and  by  increasing  their  scattering  rate.  At  high  tem¬ 
peratures  all  thermal  modes  are  excited,  but  below  the  Debye  temperature 
Td  higher  energy  modes  begin  to  freeze  out.  For  a  cubic  monatomic  crystal 
with  atomic  concentration  nA  we  have: 

Td  =  (67 v2riA) 1/3  .  (C-26) 

Because  atomic  concentrations  for  most  materials  are  comparable,  differences 
in  Debye  temperature  arise  from  differences  in  the  speed  of  sound  cs.  Typical 
values  of  the  Debye  temperature  for  thermoelectric  materials  range  from 
below  room  temperature  ( TD  =  120  K  for  Bi )  to  above  (TD  =  645  K  for  Si). 
A  hypothetical  low  temperature  thermoelectric  material  for  coolers  might 
be  able  to  gain  in  performance  by  operating  well  below  TD.  Unfortunately, 
the  small  bandgap  materials  which  are  the  best  candidates  tend  to  be  soft 
materials  with  low  Debye  temperatures,  like  Bismuth. 

Thermoelectric  performance  can  also  be  improved  by  introducing  mech¬ 
anisms  for  phonon  scattering  where  possible.  Skudderudite  materials,  which 


62 


are  good  high  temperature  thermoelectrics,  have  a  large  and  complex  unit 
cell  with  relatively  large  open  spaces.  Phonon  scattering  can  be  increased 
by  ‘filling’  these  spaces  with  otherwise  inert  atoms.  The  added  atoms  move 
relatively  freely  and  act  as  scatterers  for  phonons  in  the  basic  lattice  (see 
Mahan  1997). 
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